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Influence of disorder-induced Anderson localization and of electron-electron interaction on super¬ 
conductivity in two-dimensional systems is explored. We determine the superconducting transition 
temperature T^, the temperature dependence of the resistivity, the phase diagram, as well as the 
magnetoresistance. The analysis is based on the renormalization group (RG) for a nonlinear sigma 
model. Derived RG equations are valid to the lowest order in disorder but for arbitrary electron- 
electron interaction strength in particle-hole and Cooper channels. Systems with preserved and bro¬ 
ken spin-rotational symmetry are considered, both with short-range and with long-range (Coulomb) 
interaction. In the cases of short-range interaction, we identify parameter regions where the super¬ 
conductivity is enhanced by localization effects. Our RG analysis indicates that the superconductor- 
insulator transition is controlled by a fixed point with a resistivity Rc of the order of the quantum 
resistance Rq = h/4e^. When a transverse magnetic field is applied, we find a strong nonmonotonous 
magnetoresistance for temperatures below Tc. 
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I. INTRODUCTION 

Superconductivity |1I2) and Anderson localization [3] 
are among most important and fundamental quantum 
phenomena in condensed matter physics. These two phe¬ 
nomena are in a sense antagonists: in the case of super¬ 
conductivity the Cooper interaction creates a collective 
state with vanishing resistivity, while the Anderson lo¬ 
calization resulting from disorder-induced quantum in¬ 
terference drives the system into a state with zero con¬ 
ductivity. Therefore, when both interaction and disorder 
are present, a competition between the superconductiv¬ 
ity and localization naturally arises. This competition is 
of particular interest in two-dimensional (2D) geometry, 
where even a weak disorder makes the system an An¬ 
derson insulator. Thus, a 2D system may be expected 
to undergo a direct quantum phase transition (QPT) 
between the insulating and superconducting states, the 
superconductor-insulator transition (SIT). 

Experimentally SIT has been studied in a variety of 
2D structures, including amorphous Bi and Pb m, 
MoC [5], MoGe [7], Ta [S], InO |9lllj . NbN [T^ and 
TiN films mm, see also the reviews m- In recent 
years, there have been also a growing experimental activ¬ 
ity on SIT in novel 2D materials and nanostructures, such 
as LaAlOa/SrTiOa interfaces |16I17) . SrTiOa surfaces 
|18I19| . M 0 S 2 flakes |20l21j . FeSe thin films [22], LaSr- 
CuO surfaces [23|, and Lia;ZrNGl layered materials [24]. 
Characteristic for many of the novel structures is a strong 
screening of the Coulomb interaction due to a large di¬ 
electric constant of the substrate (such as SrTiOa). In 
addition, strong spin-orbit coupling is present in many of 
the novel materials (M 0 S 2 , LaAlOa/SrTiOa, SrTiOa). 

To drive the system through SIT, one changes a pa¬ 


rameter (film thickness, gate voltage, doping) control¬ 
ling the high-temperature sheet resistivity. With lower¬ 
ing temperature, systems with lower resistivity become 
superconducting (resistivity drops to zero), while those 
with higher resistivity get insulating (resistivity becomes 
exponentially large). The most salient observations com¬ 
mon to the majority of the above experiments are as 
follows: 

(i) Most of the experiments are interpreted as support¬ 
ing a direct transition between the superconducting 
and insulating phases, although some of them sug¬ 
gest a possibility of existence of an intermediate 
metallic phase. The critical resistivity Rc (the low- 
temperature limit of the separatrix curve separat¬ 
ing the temperature dependence of resistivity in the 
insulating and superconducting phases) is of the or¬ 
der of the quantum resistance Rq = hjAe^ ~ 6.5kD. 
However, the precise value of Rc varies from one ex¬ 
periment to another, roughly in the range between 
Rq/2 and 3Rq. 

(ii) For those systems that are superconducting (at 
low temperature T and magnetic field H), a non¬ 
monotonous dependence of resistivity on T and H 
is observed. In particular, a giant non-monotonous 
magnetoresistance is found in such systems at very 
low temperatures, T ^ T^. 

(iii) The temperature dependence of resistivity on the 
insulating side is very fast (activation or even 
stronger). 

Theoretical investigation of the interplay of interac¬ 
tion and disorder in systems with Cooper attraction has 
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a long history. Soon after the development of the micro¬ 
scopic theory of superconductivity by Bardeen, Cooper, 
and Schrieffer (BCS) [5], the question of influence of dis¬ 
order on superconductivity attracted a great deal of at¬ 
tention. It was found [25126) that the diffusive motion 
of electrons does not affect essentially the temperature 
Tc of superconducting transition, i.e the mean free path 
does not enter the expression for Tc. This statement is 
conventionally called “Anderson theorem”. 

Effects of disorder-induced Anderson localization [3] 
on superconductivity were considered in Refs. [27128) . It 
was found that, within the BCS approach, the supercon¬ 
ductivity in a disordered system persists up to the local¬ 
ization threshold and even in the localized regime near 
the Anderson transition. Furthermore, Refs. [37135] came 
to the conclusion that the mean-field transition temper¬ 
ature Tc in these regimes remains unaffected by disorder 
(i.e. the Anderson theorem holds). In a parallel line of 
research, it was discovered [33Ii3T] that an interplay of 
long-range (1/r) Coulomb interaction and disorder leads 
to suppression of T^. These ideas were put on the solid 
basis by Finkelstein [33] who developed the nonlinear 
sigma model (NLSM) renormalization-group (RG) for¬ 
malism. 

An alternative approach to the SIT known as “bosonic 
mechanism” was proposed in Refs. [33]. It takes into 
account the superconducting phase fluctuations and dis¬ 
cards completely all other degrees of freedom, in particu¬ 
lar, the localization effects. It was also proposed that an 
intermediate “Bose metal” phase may separate the su¬ 
perconductor and insulator [34) . A relation between the 
bosonic and fermionic mechanisms as a well as a status 
of the Bose metal conjecture remain quite obscure. 

Recently, Feigelman et al. [35136) found that the eigen¬ 
function multifractality near the localization threshold in 
three dimensions strongly affects properties of a super¬ 
conductor. Their remarkable finding is that Tc is dramat¬ 
ically enhanced: its dependence on the coupling constant 
is no longer exponential (as in the conventional BCS so¬ 
lution) but rather of a power-law type. This result was 
obtained on the basis of the BCS-type self-consistency 
equation, with Cooper attraction being the only interac¬ 
tion included. 

In a preceding work by the present authors m the 
influence of disorder-induced Anderson localization on 
the temperature of superconducting transition Tc was 
studied within the field-theoretical framework. Electron- 
electron interaction in particle-hole and Cooper channels 
was taken into account. The focus was put on the case of 
a weak short-range interaction (which is relevant to ma¬ 
terials with large dielectric constant, as well to cold atom 
systems). Two-dimensional systems in the weak localiza¬ 
tion and antilocalization regime, as well as systems near 
mobility edge were investigated. A systematic analytical 
approach to the problem was developed in the framework 
of the interacting NLSM and its RG treatment. The ap¬ 
proach took into account the mutual renormalization of 
disorder and all interaction constants (that, in partic¬ 


ular, leads to mixing of different interaction channels). 
This methodology allows us to explore both the cases 
of a long-range (Coulomb) interaction previously studied 
by Finkelstein [53] and of a weak short-range interaction 
within a unified formalism. More specifically, in the case 
of short-range interactions a system of coupled RG equa¬ 
tions for the problem was derived in the lowest order in 
disorder and three interaction couplings (singlet, triplet, 
and Cooper channels). 

The analysis of RG equations for the weak short-range 
interaction showed the behavior which is exactly opposite 
to that predicted by Ref. [33] for Goulomb interaction. 
It was found that the interplay of such interactions and 
Anderson localization leads to strong enhancement of su¬ 
perconductivity in a broad range of parameters in dirty 
2D systems, as well as in three dimensional (3D) systems 
near the Anderson transition (in contrast to the suppres¬ 
sion in the Coulomb case). In the latter case (vicinity of 
the Anderson transition), the microscopic theory of Ref. 
m justified previous theoretical results obtained from 
the self-consistency equation [35136) . 

This result of Ref. [37) is of fundamental importance 
and represents an unexpected physics (enhancement of 
superconductivity by localization, which is naively its 
exact antagonist). Indeed, remarkably, the localization 
physics, responsible for the increase of resistivity and 
thus driving the system towards an insulating state, fa¬ 
vors at the same time the superconductivity. The key 
condition is a suppression of the long-range component 
of the Coulomb interaction (see also Ref. [35]). This 
opens a new way for searching novel materials exhibiting 
high-temperature superconductivity: one needs the com¬ 
bination of a large dielectric background constant and 
disorder in layered structures. 

In this paper, we extend the formalism of Ref. [37] by 
deriving the RG equations to the lowest order in disorder 
but, formally, for arbitrary interaction couplings. We use 
this framework to explore systematically the interplay 
of superconductivity, interaction, and localization in 2D 
systems, with a focus on the SIT in thin films. More 
specifically: 

(i) We evaluate the temperature dependence of the 
resistivity p(T) for given bare (high-temperature) 
couplings down to the temperature Tc at which 
the finite expectation value of the superconduct¬ 
ing order parameter emerges, or else, down to the 
temperature where the system enters the insulating 
regime. 

(ii) We use the RG equations to determine the struc¬ 
ture of the phase diagram. In particular, we iden¬ 
tify parameter regions where the superconductivity 
is enhanced by localization. Our results also in¬ 
dicate that in some cases the phase diagram may 
include a critical-metal phase. 

(iii) We study the magnetoresistance near the SIT 
within two-step RG approach. Since the mag¬ 
netic held suppresses both superconductivity and 
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localization, a non-monotonous magnetoresistance 
arises, as observed experimentally. Furthermore, 
this magnetoresistance becomes very strong at low 
temperatures, again in agreement with experi¬ 
ments. Both orbital and Zeeman effects of the 
magnetic field are incorporated in the unifying RG 
scheme. 


All the above analysis is performed for the cases of short- 
ranged and long-ranged Coulomb interaction, both with 
and without spin-orbit interaction. 


The structure of the article is as follows. In Sec. El we 
introduce the NLSM formalism. The corresponding RG 
equations (valid to the lowest order in disorder and for 


arbitrary interaction strength) are presented in Sec. Ill 


The RG equations are used in Sec. to analyze the 
phase diagram in zero magnetic field. The temperature 
dependence of resistivity in zero magnetic field is dis¬ 
cussed in Sec. In Sec. VI this analysis is extended 
to calculate the magnetoresistance in a transverse and in 
a parallel magnetic field. Section m contains a discus¬ 
sion of obtained results, their implications, limitations, 
possible extensions, comparison with numerical and ex¬ 
perimental results. Finally, our results and conclusions 
are summarized in Sec. |VIII| Several Appendices contain 
technical details of the derivation of RG equations and 
of their analysis. 


II. FORMALISM 

A. NLSM action 


The action of the NLSM is given as a sum of the non¬ 
interacting part, Sa, and contributions arising from the 
interactions in the particle-hole singlet, particle- 

hole triplet, and particle-particle (Cooper), 

channels (see Refs. [321^ for review): 

s = s, + sl^^ + sl:^ + sl:l ( 1 ) 


where 




32 
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Here g is the total Drude conductivity (in units jh 
and including spin), tr = {tri 7 G 2 ,G 3 }, and we use the 


following matrices 

KL = sgnn(5„r„5“^too, 

VZm=nSnm5°‘'^tQ0, ( 2 ) 

with a,f3 = l,...,A^r standing for replica indices and 
n, m corresponding to the Matsubara fermionic energies 
En = 7rT(2n -I- 1). The sixteen matrices, 

trj — Tr Z) Sj j j — 0, 1, 2, 3, (3) 

operate in the particle-hole (subscript r) and spin (sub- 
srcipt j) spaces with the corresponding Pauli matrices 
denoted by 



Matrices tq and sq stand for the 2x2 unit matrices. 
The matrix field Q{r) (as well as the trace Tr) acts in 
the replica, Matsubara, spin, and particle-hole spaces. It 
obeys the following constraints: 

g2 = i, TrQ = 0, Q^ = C^Q^C. (6) 

The charge conjugation matrix C = it 12 satisfies the fol¬ 
lowing relation = —C. Matrix Q can be parameter¬ 
ized as g = T~^AT where the matrices T obey (symbol 
* denotes the complex conjugation) 

CT*=TC, {T-^)*C = CT-\ (7) 

In order to avoid notational confusion, it is instructive 
to compare our notation with that of the reviews |32] 
and [33] ■ In both references, a different definition of 
Pauli matrices in the particle-hole space has been used, 
namely, ir, instead of r, for j = 1,2,3. In Ref. [32] 
Pauli matrices in the spin space coincide with our def¬ 
inition (§. In Ref. [33] the spin-space Pauli matrices 
—isj (for j = 1,2,3) were used instead of our definition 
([^. The interaction terms and coincide 

with terms in Eqs. (3.9a), (3.9b), and (3.9b) of Ref. (33] 
provided the following relations between the couplings 
Fs, Ft and F^ in and and Z, Fa and 

Fc in Ref. [33] hold: Fs = —{ttv j4:)Z, Ft = {'nv/4)T2, 
and Fc = (rv/AyTc. Here the thermodynamic density of 
states V includes the spin-degeneracy factor. Note that 
Ref. |33| focuses on the case of unscreened (long-ranged) 
Coulomb interaction. Hence the interaction amplitude F^ 
in the singlet particle-hole channel is expressed through 
the frequency renormalization factor Z there. We con¬ 
sider both long-ranged (Coulomb) and short-ranged in¬ 
teractions. In the latter case the quantities Fg and Z^^ are 
independent variables. The interaction terms s[^, 
and coincide with the terms in Eqs. (3.92d), (3.92e), 
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and (3.92f) of Ref. [31] provided Fg = F^ = 
and Fc = /2. The parameters g and in are re¬ 

lated to the corresponding quantities D,Z,v of Ref. [31] 
as g = 4ttvD and Zi^ = ('ki'/4)Z and to the parameters 
G and H in Ref. [31] as g = 16/G and Z^^ = H/2 . 


B. Interaction in the Cooper channel 

The Cooper-channel interaction term can be rewritten 
as 

E E [drTr[U,LZQ] Tr[G,L“Q]. 

a,n r—1,2 j—0 

( 8 ) 


Here the matrix L“ is defined as 

(in)4=W.n<5“^<5“^^00. (9) 

However, for j = 1, 2,3 we find 

Tr [G,L“Q] =-Tr[CtJ^CLZQ] = -Tr[G,L“Q] = 0. 

( 10 ) 

(c) 

Therefore, the term describing the interaction in the 
Cooper channel is fully determined by the Cooper-singlet 
channel: 

41 = - E E / Tr [troKQ] Tr [UoKQ] ■ 

Oi,n r—1,2 

( 11 ) 


C. Relation with the BCS hamiltonian 

In general, bare values of the interaction parameters 
Fs,Ft and Fc can be estimated for a given electron- 
electron interaction U{r — r') in a microscopic hamil¬ 
tonian. It is convenient to introduce the dimensionless 
parameters 7s,t,c = ^s,t,clZui- Then their bare values 
can be written as 

Fg Ft 

7sO = -:rT4’ lto = - , , , 7c0 =-^c, (12) 

where Fg = yU{q) + Ft, 

Ft = -^^Hscr(2fcFsin(6»/2))^^^, 

Pc = Y- ^(ciscr(2A:FCOs(0/2)))^^. (13) 

Here Uscr{q) stands for the statically screened interaction 
and {■ ■ ■)FS denotes averaging over the Fermi surface. In 
the BCS case (for example, for a weak short-range attrac¬ 
tion mediated by phonons), the interaction can be writ¬ 
ten as U{r) = —{\/v)5{r) where 0 < A ^ 1. Neglecting 
screening in this case we find 

Fg « -A/2, Ft « A/2, F, « A/2. (14) 


Thus, for the BCS case (i.e. when neither screened nor 
unscreened Coulomb repulsion is taken into account), we 
get the following interaction parameters at the ultraviolet 
scale (which is given by Debye frequency cuf in the case 
of phonon-induced superconductivity): 


- 7so « 7io ~ 7co ~ -A/2. 


(15) 


If disorder is strong, ojdt 1, the relations (15) deter¬ 
mine initial values of the interaction parameters for the 
action §• In what follows, we will refer to the line de¬ 
termined by relations — 7 s = 7* = 7c as the “BCS line”. 
When disorder is weak, ojdt ^ 1, the relations (15) hold 
at the scale corresponding to the Debye frequency uif- 
Then the Cooper interaction constant is renormalized at 
ballistic scales (between ujd and l/r) such that 


-7so = 7*0 = -A/2, 

__ A/2 _ 1 

1 — (A/2) In wdt InT^'^^r 

where = ujd exp(-2/A). 


(16) 


D. F algebra and T invariance 


The NLSM action 0 involves the matrices which 
are formally defined in the infinite Matsubara frequency 
space. To perform calculations with these matrices, it is 
convenient to introduce an ultraviolet cutoff for the 
Matsubara frequencies. In addition, it is useful to intro¬ 
duce another cutoff Nm < N’j^ indicating the size of a 
non-trivial part of the Q matrix (beyond which the Q 
matrix equals A). At the end of calculations both cutoffs 
should be sent to infinity. 

Global rotations of the Q matrix with any matrix of 
the type exp(ix), where y = Ea,™ Xn^n^oo, play an im¬ 
portant role [40141] . In the limit —)■ oo and 

—)■ 0, the set of rules known as T algebra [33] 
allows one to establish the following relations (for r = 0, 3 
and j = 0,1, 2, 3): 


Tr/“Gje*^Qe-*^ 


Tr 




Tr 

8inX-ndroSjo, 


Tr TyQ -b ^ my" Tr /“tooQ 

a,n 

4 51’^^x“X-™, 


Tr 


i:uoQ 


(17) 


Using Eqs. (171, one can check that, provided Fg = —Z^, 
the action (Q is invariant under global rotations of the 
matrix Q with the matrix exp(*x) (so called F invari¬ 
ance). The constraint Tg = —Z^ corresponds to the case 
of Coulomb interaction |32|. Since the relation Fg = —Z^ 
is dictated by the symmetry of the action 0 it should 
remain fulfilled under the RG flow. 















5 


III. ONE-LOOP RENORMALIZATION-GROUP 
EQUATIONS 

A. Preserved spin-rotational symmetry 

To derive RG equations in the one-loop approximation 
(i.e., to the lowest order in disorder strength), we employ 
the background-field method and apply it to renormal¬ 
ization of the NLSM action . Details of the derivation 
can be found in Appendix In d = 2 dimensions the 
one-loop RG equations read [t = 2/['Kg)\: 


Somewhat counter-intuitively, Eq. (18 1 suggests an 


^ _ 2 
dy 


l + /(7s) + 3/(7t)-7c 


(18) 


^ =-^(1+ 7s)(7s+ 37* + 27c+ 47 ^), (19) 

^ = -^(1 + 7t) Is-It- 27c( 1 + 27t - 27c) 


( 20 ) 


^ = -27c - h [(1 + 7c)(7s - 37t) - 27 ^ -b 47 ^ 
dy 2 L 


+ 67 c(| 7 t-ln(l-b 7 t)) , (21) 


d In Zcj t 

dy 2 


(is + ^It + + 47c) , 


( 22 ) 


where y = \n{L/l) {I denotes the mean free path) and 
f{x) = 1 — (1 -b 1/a;) ln(l -b a:). These RG equations de¬ 
scribe the evolution of the system with spin-rotational 
and time-reversal symmetries upon changing the char¬ 
acteristic length scale L. We stress that RG equations 
(181 - ( [ 2 ^ satisfy the particle number conservation since 
dJ^ui+^s)/dy = 0. Further, it is worth emphasising that 
the right-hand-sides of the equations are nonsingular in 
the limit of Goulomb interaction, 7 ^ = —1. 

The ultraviolet value of the NLSM coupling t that de¬ 
scribes the disorder strength is given by the dimension¬ 
less Drude resistivity. The renormalization of t at larger 
scales involves the contributions to the resistivity induced 
by interference effects and by virtual (elastic) processes 
due to interactions in particle-hole singlet ( 7 ^) and triplet 
( 7 t), as well as in Cooper channel ( 7 c). 

We emphasize that Eqs. (18 1 - ([22|) are obtained in the 


lowest order in t but they are formally exact in interac¬ 
tions worth noting that the Cooper-interaction 

coupling 7 c enters all the RG equations only in a polyno¬ 
mial way. Interestingly, the contribution of Cooper chan¬ 
nel to the renormalization of t is fully described by the 
linear term only, thus rendering Eq. (18) for arbitrary 7 c 
the same as in the weak-coupling limit 52]; l 7 c| ^ 1 - 


The first term in Eq. (21) describes the standard BCS 


instability; in accordance with the “Anderson theorem” 
this term is not affected by disorder. Moreover, the “An¬ 
derson theorem” manifests itself in Eq. (21) through the 


absence of the terms tO^jc) on the right hand side. To 
the lowest order in interaction couplings, the effect of 
disorder on the renormalization of 7 c is solely due to the 
presence of the interaction in the particle-hole channels. 


insulating behavior (an increase of the resistivity with 
increasing L) for 7 c —>■ — 00 . We note, however, that 
the (dimensionless) physical resistivity p is not exactly 
equal to the NLSM coupling t because of the inelastic 
contribution to the conductivity governed by supercon¬ 
ducting fluctuations, see Sec. |V] below for details. Near 
the superconducting instability (for large [yd ^ 1 ), this 
antilocalizing inelastic contribution to the conductivity 
becomes large. 

Furthermore, towards the superconducting instability, 
7 c —> — 00 , the disorder-induced renormalization of 7 c in 
Eq. (21) is dominated by the term —2tj^ which tends 


to impede a development of the superconducting insta¬ 
bility. Thus, if Eqs. (18) - (22) would constitute the ul¬ 


timate truth, the superconducting instability would not, 
strictly speaking, develop. An explanation for this ap¬ 
parent paradox is as follows. It turns out that the one- 
loop RG equations become insufficient in a vicinity of 
the superconducting instability, namely, on scales larger 
than Lx where [yd reaches a value ~ 1/t ^ 1. In other 
words, the weak-disorder condition of validity of the one- 
loop RG, t 1, should in fact be supplemented by the 
condition flyd ^ 1 . 

The emergence of the latter condition (and thus of 
the scale Lx) becomes evident from a comparison of the 
terms of the zeroth and the first order in t in Eq. (21 1 . 


This scale Lx arises also in the calculation of the conduc¬ 
tivity (see Sec. 0: at this scale the inelastic contribution 
to the conductivity reaches in magnitude the elastic one. 
We expect that in the vicinity of the superconducting 
instability higher-loop terms of the type t{t')c)^ in the 
beta-function for t and equation govern¬ 

ing renormalization of 7 c should emerge. Upon resumma¬ 
tion, they are expected to restore the divergence of yc at a 
scale Lc slightly larger than Lx- At the same time, since 
the second-loop (fc = 2 ) terms are similar to those de¬ 
scribing mesoscopic fluctuations of the superconducting 
order parameter |55I5H| . we expect for jydt > 1 (i.e., for 
temperatures slightly above the transition) strong spa¬ 
tial fluctuations of observables (in particular, of the local 
tunneling density of states |45I46) . as observed in exper¬ 
iments, see, e.g., Ref. my 

To the lowest order in yc, Eqs. (18) - (22) coincide 
with results obtained by Finkelstein long ago [52|. Re¬ 
cently, one-loop RG equations beyond the lowest order 
in interactions were reported in Ref. [53] for the case of 
preserved spin-rotational and time-reversal symmetries. 
It should be stressed, however, that our RG equations 
- ( [ 2 ^ differ from those of Ref. [53] . It is instructive 
to highlight the difference. First of all, the right hand 
side of a RG equation for y^ in Ref. [53] [see Eq. (A12) 
there] contains a term proportional to ty/ rather than 
to t{l + 7s)7c as in our Eq. (19). Since the quantity 
Zuj + Ts = Zi^(l -|- 7 s) should have no renormalization 
by virtue of the particle number conservation, this would 
imply the presence of a term proportional to t7c/(l +7s) 
in the RG equation for Z^^. Being divergent for the case 
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of Coulomb interaction, 7 s = —1, such a term would, 
however, violate the J^-invariance of the NLSM action 
0 and is thus not allowed. Second, the RG equation 
for 7 t reported in Ref. [13] does not contain the term 


proportional to ^ 7 ^, in contrast to our Eq. (|20|). Finally, 
the RG equation for 7 c reported in Ref. 


43 


contains 

an additional term proportional to lii(l + 7s) as com¬ 
pared to our Eq. (21). We note that a similar term was 
reported by Belitz and Kirkpatrick in Ref. [33] [see Eq. 
( 6 . 8 g) there]. In our opinion, such terms, divergent for 
the case of Goulomb interaction, 7 s = —1, cannot appear 
in the course of renormalization of J^-invariant operators, 
including Sc- In Ref. m, the appearance of a term pro¬ 
portional to t 7 cln(l -|- 7 s) in the RG equation for 7 c of 
Ref. |39| was attributed to an improper treatment of the 
gauge invariance. In our background-field RG calcula¬ 
tions, terms proportional to ln(l -|- 7 s) do appear in the 
course of renormalization of Ec at intermediate steps but 
cancel each other in the final results, in agreement with 
the J^-invariance, see Appendix [X] 


B. General case 


The RG equations (18) - (22) have been derived for the 
case of preserved spin-rotational symmetry. We are now 
going to generalize them to systems with spin-rotational 
symmetry broken (partly or fully) due to spin-orbit cou¬ 
pling and/or spin-orbit impurity scattering. Both these 
symmetry-breaking mechanisms induce finite relaxation 
rates (I/rf, 1/t^, I/t"/) for corresponding components 
of the electron spin. The relaxation rates determine 
the mass of the corresponding triplet modes (diffusons 
and cooperons). As an example, the mode correspond¬ 
ing to the spin component acquires a mass propor¬ 
tional to 1/rf -I- I/r/. This mode thus become effectively 
frozen and drops out of RG equations at length scales 
L » if ~ [l/iDry) + l/(I4r/)]-i/2. 


In the presence of spin-orbit coupling, there is the spin 
relaxation due to D’yakonov-Perel’ mechanism. The cor¬ 
responding relaxation rates are given by 1 /r/’*'’^ ~ 
where Ago denotes the spin-orbit splitting [48]. There¬ 
fore, all triplet modes (both for diffusons and cooperons) 
are suppressed at the length scales L ^ Lgo = vp/^so, 
i.e., the number of triplet modes contributing to the RG 
equations is n = 0. In the case of a 2D electron sys¬ 
tem with the spin-orbit impurity scattering but with¬ 
out spin-orbit coupling, the spin relaxation is anisotropic: 
1 /r/ = I/tso, = 0 , where I/tso denotes the skew 

scattering rate [49]. Thus, for L Lg = the 

triplet modes corresponding to the total spin component 
Sz remain massless. Therefore, in this case n = 1 triplet 
mode still contributes to the RG equations. 

If the spin-orbit coupling and spin-orbit scattering are 
both present, then different regimes with n = 3, n = I 
and n = 0 can be realized depending on the relations 
between L, Lgo and Lg. For all three cases the one-loop 


RG equations can be written as 


^ _ 2 
dy 


n — 1 


+ fils) + nf^jt) - 7c 


(23) 


^ = -^(l + 7s)(7s+’^7t + 27c + 47/), (24) 

^ = -^(1 + 7t) ls-{n- 2)7t 


- 27 c (1 -I- 27 * - 27 c) 


(25) 


^ = -27c “ ^ “ 27c + 47c 


din Ziz 
dy 


+ 2n'yc{lt - ln(I -I- 7 t)^ 
= ^ (7s + n'yt + 27 c -f 47/) . 


(26) 

(27) 


In the case n = 0, Eq. (25) should be omitted. The RG 


equations (23) - (26) constitute one of the main results of 


the paper. In the rest of the paper, we will analyze these 
equations to investigate phase diagrams and observables 
for the cases of preserved and broken spin-rotational sym¬ 
metry. 


The system of RG equations (231 - (26) has the fixed 


plane 7 ^ = — 1 corresponding to the case of long-ranged 
Goulomb interaction. In fact, this statement is not re¬ 
stricted to the one-loop RG equations. The existence of 
such a fixed plane is a consequence of the particle-number 
conservation and of the J^-invariance of the NLSM action 
0 . Due to the charge conservation, RG equations for 7 ^ 
and z are related to all orders in t: 


dy 


■ 7s)Cz) 


dZgj _ rz c 
—1 
dy 


(28) 


The value of the anomalous dimension C,z at a fixed 
point determines the dynamical critical exponent z = 
d + (*. The latter controls the temperature behavior of 
the specific heat, Cy ^ [SD]. Typically, one expects 

that z ^ d {Q < 0 ) which implies the instability of the 
fixed point in the plane 7 ^ = — 1 with respect to the 
increase of 75 . 

It is worth reminding the reader that RG equations 


(23) - (27) are of one-loop order with respect to diffusive 
modes (i.e., are derived by expansion of the right hand 
side to the lowest nontrivial order in t) but are exact 
in interaction. Typically, one expects that one-loop RG 
equations are valid until entering the insulating (strong- 
disorder) phase, i.e. for t < 1. This requires a tacit 
assumption that in the expansion of the right hand side 
of RG equations in powers of t all coefficients (which are 
functions of interaction amplitudes) are of the order of 
unity. In the case of superconducting instability, 7 c di¬ 
verges at some scale Lc, so that coefficients of the expan¬ 
sion in powers of t become much larger than unity. As 
discussed in Sec. [Ill A[ near the superconducting insta¬ 
bility (i.e., at | 7 c| ^ I) the general condition of validity 
of the one-loop approximation t I becomes more re¬ 
strictive: t| 7 c| < 1. Similarly, near the Stoner instability 
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(which corresponds to the divergence of 74 ) the two-loop 
analysis [51152] demonstrates that expansion in t is jus¬ 
tified for t < 1/74 1 . 

Up to now we have discussed the renormalization as a 
flow of couplings with the length scale. In practice, one 
usually has a sufficiently large system and the infrared 
cutoff is controlled not by the system size but rather by 
the temperature T. In this situation, the renormaliza¬ 


tion due to the contributions to RG equations (23) - (27) 


induced by interactions should be stopped at the lengt h 
scale L-t which is determined as follows [see Eq. (A15)]: 


rp _ \ ( 

T \Lt ) t{LT)Z^{LT)' 


(29) 


where to = ^(0 — Zuiil)- This transformation 

of temperature into the length scale [SSI allows us to in¬ 
vestigate the temperature dependence of observables. In 
particular, the electrical resistivity in the absence of mag¬ 
netic field is addressed in Sec. jV] The inclusion of mag¬ 
netic field induces two additional length scales, Ih and 
Iz, related to the orbital and Zeeman effect of magnetic 
field and leading to the magnetoresistivity. Sec. lYg 


IV. PHASE DIAGRAM AT ZERO MAGNETIC 
FIELD 

A. Preserved spin rotational symmetry 


We start our analysis of RG equations (23) - (27) from 


the case in which spin rotational and time reversal sym¬ 
metries are preserved, i.e., there are n = 3 triplet modes. 
We note that in notations of Ref. [35] this case is termed 
as G(LR) for Coulomb interaction and G(SR) for short- 
ranged interaction. 


1. Coulomb interaction 

For the case of Coulomb interaction, js = which 
is the fixed plane of Eqs. (23) - (26), the RG equations 
can be simplified as (we set n = 3) 


* _ 2 

dy 

t 


2 + 3/(74) - 7 c 


(30) 


(1 + 74) (1 + 7i + 27 c( 1 + 274 - 27c)) , ( 31 ) 


djt 
dy 

^ = -272 + ^ (1 + 7 c)(l + 374) + 27 c '(1 - 27 c) 
dy 2 L 


t r 


670 ( 7 * - ln(l -f 74 )) 


(32) 


Let us now analyze fixed points of Eqs. (30) - (32). It 


turns out that the structure of the set of fixed points 
and of the three-dimensional phase diagram is very rich. 
Specifically: 


• There is a marginally unstable line of fixed points at 
t = 7c = 0 (with arbitrary 74). These fixed points 
describe a conventional clean Fermi liquid without 
Cooper-channel attraction. 

• There is a line of fixed points at t = 0 and 7 c = —00 
(with arbitrary 74) corresponding to the supercon¬ 
ducting (SC) phase. 


Further, Eqs. (30) - (32) contain also the attrac¬ 


tive line of fixed points at 74 = 00 and 7 c = 1. The 
divergence of 74 corresponds to a ferromagnetic in¬ 
stability. 


Formally, in Eqs. (30) - (32), there exists also a 
fixed point at 74 = —1, 7 c = 0 and t = 00 . While 
the range of t > 1 is beyond the accuracy of the 
one-loop RG, it is expected on general grounds that 
full RG equations should contain an attractive fixed 
point (or a family of fixed points) with t = 00 de¬ 
scribing the insulating phase. 


• Within Eqs. (30) - (32) there is a possibility at 


some length scale to enter the phase with 74 = — 1 . 
At this length scale there are finite values 7 c < 0 
and t. We note that 74 = —1 corresponds to the 
infinitely strong attraction in the triplet particle- 
hole channel indicating a possibility of exciton con¬ 
densation. Since the value 74 = —1 is reached at 
a length scale close to Lx, full RG equations are 
needed to study a competition of exciton conden¬ 
sation in the spin channel and superconductivity in 
the Cooper channel. We leave this as a prospect for 
future research and do not discuss a possibility of 
exciton condensation in the rest of the paper. 


Going beyond the one-loop RG equations (|30|) - 

1 , 


we expect a fixed point at 7 c = — 00 , t 
and a certain value of 74 governing the transition 
between the superconductor and insulator phases. 
The corresponding phase boundary is a critical sur¬ 
face with a flow towards this SIT fixed point origi¬ 
nating at the trivial fixed point with t = 7 c = 7 * = 
0. We will discuss the SIT fixed point in more detail 
in Sec. IVIII below. 

• Similarly, we expect strong-coupling fixed points 
that control the ferromagnet-insulator and the 
ferromagnet-superconductor transitions. We will 
not discuss these fixed points in the present paper 


Let us now discuss properties of the emerging phases 
(see Figs. and and corresponding fixed points, in 
more detail. 

Superconducting phase. We first note that within the 
RG equations (30) - ([32|) the superconducting line of 


fixed points at t = 0 and 7 c = —00 is unstable, which 
makes the superconducting phase formally unreachable. 
As we have already discussed, this indicates a failure of 




























FIG. 1: (Color online) The case of preserved spin rotational 
symmetry with Coulomb interaction, 7 s = —1: the RG flow 
obtained from numerical solution of Eqs. (|30[) - (32l. The 


initial condition Axes 7 to = 0.2. The arrows indicate the flow 
towards the infrared. The gray region indicates the part of 
the phase diagram which is not accessible within one-loop 
RG equations. The lines describing the flow to the supercon¬ 
ducting (SC), insulating (I), and ferromagnetic (F) phases are 
shown in red, green, and blue, correspondingly. Gray flow 
lines correspond to the region of superconductor-insulator 
transition (SIT). 


the one-loop (lowest order in t) RG equations near the 
superconducting instability. In the absence of disorder 
(i.e., at t = 0), Eq. (32) describes the usual BCS-type 
scenario. The Cooper-channel interaction 7 c diverges at 
some finite length scale L,, as 7 c(T —>■ L^) ~ — l/(Lc~^)- 
To estimate the length scale Lc in the case of finite dis¬ 
order, we shall use the scale Lx defined by the condi¬ 
tion \jciLx)\ = ^/t{Lx) ^ 1. Assuming that the di¬ 


vergence of 7 c is of the BCS type, we get an estimate 
(Lc — Lx)/Lx ~ t{Lx) ^ 1- Thus, while the one-loop 
RG is not sufficient to follow the flow up to the singu¬ 
larity scale Lc, it works up to a scale Lx which is only 
slightly smaller than Lc. 

Insulating phase and superconductor-insulator transi¬ 
tion. On general grounds, we assume that once the RG 
flow reaches t ~ 1 , the system is in the insulating phase, 
i.e., it flows into the insulating (I) fixed point with t = oo. 
On the other hand, as discussed above, if t remains small 
when | 7 c| reaches a value 1/t, the system flows into a su¬ 
perconducting fixed point. There should be thus a fixed 
point at t ~ 1 (i.e., with resistivity of order of quantum 
resistance Rq) and certain values of 7 c and 7 t that con¬ 
trols the quantum phase transition between superconduc¬ 
tor and insulator, see Sec. VII for a further discussion. 


At small values of t and 7 c < 0, 7 t > 0, the separatrix 
surface between the two phases is parametrized by the 
following equation; t = 47 c/(I + ^It)- 


0.5 


N 

+ 


-0.5 



-1 


-0.5 

rco/^ 1 + 7co 


FIG. 2: (Color online) The case of preserved spin rotational 
symmetry with Coulomb interaction, 7 ^ = — 1: a projection of 
the phase diagram on the 'yto — 7 c 0 plane. The color indicates 
the value of the Drude resistivity to at which the quantum 
phase transition from SC to I occurs. Above the dashed line, 
the FM phase appears in addition to SC and I phases. The 


figure is obtained from numerical solutions of RG Eqs. (30) - 

(^. 


“Ferromagnetic” phase. For the attractive line of fixed 
points at 7 t = 00 and 7 c = 1 , the value 7 c = 1 is fixed 
by a cancelation of terms in the right-hand side of Eq. 
(32) which are proportional to 74 1. The divergence 


of 7 t occurs at some finite length scale LpM- Due to a 
delocalizing effect of the interaction (Altshuler-Aronov) 
contribution to renormalization of the resistance at large 
7 t, the fixed point value t[LpM) remains finite and is 
non-universal (i.e., determined by the initial conditions). 
Therefore, Eqs. ([30| - ([32|) predict ferromagnetic metallic 


phase with a non-universal resistivity. Strictly speaking, 
one-loop equations are insufficient to describe accurately 
the regime t'jt ^ 1 (see Refs. |51I52) ) but this is not ex¬ 
pected to modify essentially the emergence of instability. 

However, since the emergent fixed points are character¬ 
ized by a finite value of dimensionless resistivity t{LpM), 
the diffusive RG continues at larger scales. Specifically, 
to describe properly the system at scales larger than 
LpMi one needs to take into account breaking of spin ro¬ 
tational symmetry and derive a new set of RG equations. 
In this case all triplet diffusive modes in the particle- 
hole channel and singlet and triplet modes in the Gooper 
channel are suppressed. One can thus assume that the 
system at L > LpM is described by RG equations ([^ 
with n = 0 , 7 c = 0 and 7 s = —1 that results in in¬ 
sulating behavior at large length scales. Moreover, due 
to enhanced spin fluctuations near the Stoner instability, 
the system at L > LpM can demonstrate a spin-glass 
behavior [^. In what follows, we shall term this phase 
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FIG. 3: (Color online) The case of preserved spin rota¬ 
tional symmetry with Coulomb interaction, 7 ^ = —1: depen¬ 
dence of t (“renormalized Drude resistivity”) on the length 
scale across the quantum phase transition between super¬ 
conducting (SC, red curves) and insulating (I, blue curves) 
phases. The curves are obtained from numerical solutions of 
RG equations (|^ - ([32l) for 7 co = —0.25, qto = 0.01 and 
to = 0.05, 0.1, 0.12, 0.14, 0.15, 0.18, 0.2, 0.22 (from the 
bottom to the top). 


ferromagnetic (FM) for simplicity. 

Overall RG flow and phase diagram. A part of the RG 
flow for Eqs. (30) - (32) is shown in Fig. In general, a 
projection of the flow in a three-dimensional parameter 
space onto a 2D plane, as in Fig. depends on initial 
conditions for the couplings. For the plot shown in Fig. 

we have assumed a realistic relation between the triplet 
(third axis) and Cooper amplitudes, which has allowed 
us to avoid intersections in the projected flows. Further¬ 
more, the RG flow is shown only in the region of validity 
of the one-loop approximation: tmax{l, lyd} ^ 1. The 
flows towards the superconducting, insulating, and fer¬ 
romagnetic phases are plotted in red, green, and blue, 
correspondingly. The grey part of the flow describes the 
vicinity of the SIT. One of the grey curves is the separa- 
trix between the superconducting and insulator phases. 
However, the one-loop precision is insufficient to deter¬ 
mine the separatrix in the region tmax{l,| 7 c|} > 1 . 
At small values of 7 ^ the separatrix is parametrized by 

^ = 47c/(l + 37t)- 

The phase diagram expected on the basis of the RG 
equations ( [^ - (32) is shown in Fig. in the plane 
of bare interaction couplings 7 co and 7 * 0 . For 7 co < 0, 
the superconducting phase exists at small values of tg. 
For given 7^0 and ■jto the quantum phase transition from 
superconductor to insulator occurs with increase of tg- 
In addition, for a sufficiently large 7 tg (above the dashed 
line) a ferromagnetic phase emerges. In this part of the 
7 c 0 ~ 7 io plane, a sequence of transitions S - I - FM 
- I takes place with increasing bare resistivity tg. For 


FIG. 4: (Color online) The case of preserved spin rotational 
symmetry with Coulomb interaction, 7 s = —1: dependence 
of t on the length scale across the quantum phase transi¬ 
tion between superconducting (SC, red curves), insulating (I, 
green curve), and ferromagnetic (FM, blue curves) phases. 
The curves are obtained from numerical solutions of RG 


equations (301 - (32) for 7^0 = —0.1, 7 to = 0.4 and to = 
0.015, 0.0165, 0.18325, 0.18326, 0.184, 0.02, 0.022, 0.025 
(from the bottom to the top). For higher values of to, another 
insulating phase (not shown on this scale) emerges. Inset: de¬ 
pendence of and Tfm on to- 


7 cg > 0 , there is no superconducting phase; changing tg 
drives a transition from the ferromagnetic to the insulator 
phase. 

The dependence of the NLSM coupling t on the length 
scale L across the quantum phase transition from the 
superconducting to insulating phase (in the part of the 
phase diagram in Fig. where FM phase does not oc¬ 
cur) is shown in Fig. This dependence dominates the 
corresponding evolution of the total electrical resistivity 
p (apart from a narrow region close to the superconduct¬ 
ing instability, where the inelastic contributions due to 
fluctuating Cooper pairs becomes dominant, see Sec. [V] 
for details). 

In Fig. I^we choose the values of 7 cg and 7 tg such that 
the FM phase exists in addition to the SC and I ones. We 
thus show the length dependence of t across the quantum 
phase transitions from SC to I and from I to FM phases. 
We note that within RG Eqs. (30) - (|32|) the insulat¬ 


ing phase (between SC and FM phases) exists in a very 
narrow interval of tg, see Fig. As one can see, the 
scale Lx (at which red curves in Fig. are stopped), 
which yields approximately the superconducting coher¬ 
ence length, is larger than the BCS coherence length 
= Z exp(—l/ 27 cg). In the ferromagnetic phase, the 
corresponding length scale Lpm (where blue curves end) 
is still larger than Lx- 

At finite temperature, the interaction contributions to 


the RG equations (30) - (32) are stopped at the length 
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scale Lt- Neglecting the difference between Lt and the 
temperature-induced dephasing length (which cuts off 
the localization corrections), we can stop the whole RG at 
Lt- Then the transition temperatures to superconduct¬ 
ing (Tc) and ferromagnetic phases {Tfm) is esti mated 
as follows (see also a discussion in the end of Sec. fill: 
Tc « (l/T){l/Lxf and TpM « LpAif ■ A typi¬ 

cal dependence of Tc and Tfm on to is shown in the insets 
to Figs. |^and|^ The effect of disorder on Tc depends on 
the sign of the term in the square brackets in the right 
hand side of Eq. (32). It occurs that for 7 c < 0 and 
7 t > — 1 this term is always positive, except for a small 
region at small negative values of 7 c and —1 <74 < —1/3. 
Therefore, as was first found by Finkelstein [32], disor¬ 
der in the presence of Coulomb interaction suppresses 
the superconducting phase (i.e., lowers Tc). At the same 
time, disorder induces the ferromagnetic phase which ex¬ 
ists in an intermediate range of disorder. This implies a 
nonmonotonous dependence of Tfm on to- 

We note that Tc evaluated from the RG equations (30) 
- (321 is in fact somewhat larger than the true supercon¬ 
ducting (Berezinskii-Kosterlitz-Thouless) transition tem¬ 
perature Tbkt due to the presence of phase fluctuations 
of the order parameter at temperatures below Tc, see Sec. 
I^for more detail. The relative difference between Tc and 
Tbkt is, however, small for weak disorder, and thus does 
not essentially affect a much stronger variation of Tc with 
disorder explored in this paper. 


2. Short-ranged interaction 


In the case of short-ranged interaction, RG equations 


(23) - (26) with n = 3 read 


dt 

dy 

d^s 

dy 

d'jt 

dy 

d'ye 

dy 


1 +/(7s)+ 3/(74)-7c 


- 2 (1 + 0''*) + ^0'* + 27c + 47c): 


(33) 

(34) 


2(1 + 7*) -7s+7* + 27c(1 + 274 - 27c) , (35) 
-27c + ^ (1 + 7c)(-7s + 374) -k 27^(1 - 27c) 


67c (7* - ln(l +7t) 


(36) 


Contrary to the Coulomb-interaction case (where we 
had 7 s = — 1 ), the singlet particle-hole amplitude 7 s is 
not fixed now, so that the RG flow occurs in the four¬ 
dimensional parameter space. However, the structure 
of the set of attractive fixed points (quantum phases) 
and of fixed points describing quantum phase transitions 
between them remains qualitatively the same as in the 
Coulomb case. Specifically, the fixed points of the RG 
flow for the short-ranged interaction are as follows [55] : 


• There is a surface of clean-Fermi-liquid fixed points 
at t = 7 c = 0 (with arbitrary 74 and 7 s). 


• The fixed-point surface at t = 0 and 7 c = —00 
corresponds to the superconducting phase. 


• The line of fixed points with 7 s = — 1, 74 = 00 , 7 c = 
1 , and arbitrary t, is attractive in the js direction. 
Therefore, the RG equations (33)-(36) lead to the 
same ferromagnetic phase that exists in the case of 
Coulomb interaction. 


• Exactly as in the Coulomb case, there should be a 
fixed point (or a family of fixed points) with t = 00 
describing the insulating phase. 

• Eor the same token as in the Coulomb case, a SIT 
fixed point with t ^ 1 should separate the super¬ 
conducting and insulating phases. 


The phase diagram for a given 7 so > — 1 is similar 
to that for the case of Coulomb interaction, 7^9 = — I 
(shown in Fig. [^. With increase of 7 ^ 0 , the destruction 
of the superconducting phase gets shifted towards larger 
values of to- The crucial difference between the cases of 
short-ranged and Coulomb interactions is the existence 
of large region of the phase diagram with Lx < 

(and thus Tc > In the case of a bare repulsion 

in the particle-hole channel, 7^0 < 0 and 740 > 0 , the 
superconducting transition temperature is typically lower 
than the clean BCS result, T^ < T^^^ (see Fig. [^. 
However, the situation changes if the bare interaction in 
the triplet particle-hole channel is attractive, 740 < 0 . 
As illustrated in Fig. a significant part of the phase 
diagram is occupied by superconductor with T^ > 

It should be emphasized that the superconducting phase 
with enhanced Tc exists also for 740 > 0. However, it 
occurs only in a small region of 740 , | 7 co|, l 7 so| 1 (see 
Fig. §. Typical RG evolution of the resistance t in 
this region of initial values of interactions is shown in 
Fig. i Being initially suppressed by disorder, Tc can 
be significantly (several orders of magnitude) enhanced 
with respect to Tj^^^ near the superconductor-insulator 
quantum phase transition, as illustrated in the inset to 
Fig. This is in agreement with the conclusion of our 
work [37] where RG equations (18l-(21) with the right- 
hand sides expanded to the lowest nontrivial order in 7 s, 
74 and 7 c were analyzed. 

The mechanism of enhancement of the transition tem¬ 
perature is as follows. For small initial values of inter¬ 
action parameters | 7 so|, l 7 to|, l 7 co| ^ to 1 , the renor¬ 
malization of the Cooper interaction amplitude occurs in 
two distinct steps. At the first step of the RG flow, the 
interaction is renormalized due to the presence of disor¬ 
der (the terms proportional to t), while at the second 
step the standard BCS-type renormalization (the term 
— 27 /) takes place. At the first step of renormalization, 
we can linearize the RG equations (33)-(36) in interac¬ 
tion parameters and neglect the term — 27 /. Then, in the 
course of RG flow, the interaction amplitudes approach 
the BCS line 7s = 74 = —7c, converting the repulsion 
in singlet and triplet particle-hole channels into attrac¬ 
tion. This is the consequence of the (weak) multifractal- 























11 



7cO 



FIG. 5: (Color online) The case of preserved spin rotational 
symmetry with short-ranged interaction, 7^0 = —0.05: a pro¬ 
jection of the phase diagram on the 7 to — 7 c 0 plane. The color 
indicates the ratio of Tc/{Tc + for to = 0.06. The 

dashed curve separates the regions with Tc < and with 

Tc > . The figure is obtained from numerical solutions 

of RG Eqs. ([Mil - (fel. 


ity of the noninteracting fixed point. At some scale Li 
such that hiLi/l = 1 /to ~ 1 /^ the interaction couplings 
become of the order of the resistance: | 7 s 7 ,c| t. Pro¬ 
vided ItsoUtmUTcoI ^ to, the resistance at this scale 
t{Li) ~ t^/max{|7so|, |7tol, ITcoII < 1 and all interac¬ 
tion parameters are still much smaller than unity. After 
the length scale Li the second step of RG flow starts, 
where in Eq. (21) one can neglect terms proportional 
to t compared to the disorder-independent term — 27 /. 
Thus, the Cooper interaction 7 c flows according to the 
standard BCS RG equation for a clean system with the 
initial value 7 c(Ti) ~ t(Li) rather than 7 co. Hence, 
we find the following rough estimate for the transi¬ 
tion temperature: Tc ^ (l/r)(^/Li)^ exp(— 1 /| 7 c(Ti)|) ^ 
(1/r) exp(—2/to) 3> Tc^^ (see Appendixfor details). 

In short, the role of the non-interacting disorder- 
induced multifractality is to enhance the interaction in 
the Cooper channel such that it becomes comparable 
to the resistance. After that, the divergence in the 
Cooper channel is driven by the standard mechanism 
(the same as in a clean system). The enhancement of 
Tc occurs in an intermediate range of disorder (to be¬ 
tween | 7 io| and jy^ol^^^). For a weaker disorder to < 
max{| 7 co|, | 7 to|j l 7 so|}j one can find a suppression of the 
transition temperature instead of enhancement, see solid 
curve in the inset to Fig. 

If the disorder scattering rate l/r exceeds the Debye 
frequency lud, the starting point of the RG flow will be 
likely located not fa r from the BCS line, j^o = —Jto = 
—7c0 (see Sec. IIC). For such initial conditions, the de- 


FIG. 6 : (Golor online) The case of preserved spin ro¬ 
tational symmetry with short-ranged interaction, 7 so = 
—0.05: dependence of t on the length scale across the 
quantum phase transition between superconducting (SC, red 
curves) and insulating (I, blue curves) phases. The curves 
are obtained from the numerical solution of RG equations 
(331 - (|M| for 7 co = —0.04, jto = 0.005, and to = 


0.01, 0.03, 0.05, 0.06, 0.07, 0.078, 0.085, 0.095, 0.105 (from 
bottom to top). For higher values of to, an insulating phase 
(not shown on this scale) emerges. The dependences of 
Tcl{Tc -\- on to for 7^0 = —0.04, 7*0 = 0.005 and 

7 so = —0.05 (solid curve) and for 7 co = — 7 s 0 = 7 to = —0.1 
(dashed curve) are shown in the inset. 


pendence of Tc/{Tc + on to is shown in the inset 

to Fig. [^ by dashed curve. It is worth stressing that 
in the case of initial interaction parameters on the BCS 
line there is no initial decrease of transition temperature 
with increase of to- This is because the second term in 
the right hand side of Eq. (36 1 is negative on the BCS 
line for 0 > 7 c > —0.41. The dependence of Tc/T^f'®*" on 
7 co and to on the BCS line is shown in Fig. 

Let us now turn to the region in the phase diagram. 
Fig. where the ferromagnetic (FM) phase emerges. 
Figure]^shows typical dependences of resistance t on the 
length scale L across the transition from SC to FM phases 
with increasing to- The inset presents dependences of 
both critical temperatures (Tc and Tfm) on to. 


B. Broken spin rotational symmetry 


In the presence of spin-orbit coupling and/or spin-orbit 
scattering the spin-rotational symmetry is broken. At 
length scales L ^ max{Tso, T^} the spin-rotational sym¬ 
metry is completely broken and all triplet modes are sup¬ 
pressed such that n = 0. In Ref. [39] this case is referred 
to as SO(LR) for Coulomb interaction and SO(SR) for 
short-ranged interaction m- 
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FIG. 7: (Color online) The case of preserved spin rotational 
symmetry with short-ranged interaction on the BCS line: the 
color density plot for the ratio Tc/(Tc -|- in the to — 7co 

plane. The dashed black lines separate the regions with Tc > 
and Tc < T^^^. 


1. Coulomb interaction 


For n = 0 and for the case of Coulomb interaction, 
7 s = —1, the one-loop RG equations, Eqs. (23) - (26), 
take the form 


dy V 2 

t ' 
2 . 


dy 


1-f 7c + 27^(1 - 27c) 


(37) 

(38) 


The first, perturbative study of the effect of interaction 
on conductivity of a disordered system in the presence of 
spin-orbit scattering has been performed in Ref. [53] . To 
the lowest order in 7 c Eqs. (37) - (1^ coincide with the 


one-loop RG equations derived in Refs. 

Since the spin-orbit interaction kills the contribution of 
the triplet channel, while the particle-hole singlet ampli¬ 
tude remains fixed, 7 s = — 1, the RG flow now occurs in 
a 2D parameter space, t and 7 c. The structure of phase 
diagram is governed by the following fixed points: 


Equations (37) and (38) possess a clean-Fermi- 


liquid fixed point at t = 7 c = 0 which is marginally 
unstable. 


• There is the fixed point at t = 0 and 7 c = —oo 
corresponding to superconducting (SC) phase. 

• As in all other symmetry classes, there should be 
the insulating (I) phase with t = oo. It is, however, 
not reachable within the one-loop RG equations. 


FIG. 8: (Color online) The case of preserved spin rota¬ 
tional symmetry with the short-ranged interaction, 7^0 = 
—0.05: dependence of t on the length scale across the 
quantnm phase transition between snpercondncting (SC, red 
curves) and ferromagnetic (FM, blue curves) phases. The 
curves are obtained from the numerical solution of RG equa¬ 
tions (|33l - ([so]) for 7co = —0.1, 7to = 0.2, and to = 
0.02, 0.04, 0.05, 0.063, 0.065, 0.07, 0.075, 0.08, 0.09, 0.1 


(from bottom to top). The dependences of TcjTf^^ and 
TpM iTc^^ on to are shown in the inset. 


phase. This fixed point appears at the borderline 
of applicability of one-loop RG equations, t* ^ 1, 
so that we do not have a rigorous argument in 
favour of existence of the CM phase. We find, how¬ 
ever, very plausible that the attractive character 
of this fixed point is not destroyed by going be¬ 
yond one loop. The emergence of this fixed point 
can be traced back (i) to the competition of weak 
antilocalization (enhanced by delocalizing effect of 
repulsive Cooper-channel interacti on) with the lo¬ 
calizing Coulomb repulsion in Eq. (37), and (ii) to 
the competition between Cooper instability and the 
disorder-induced suppression of the interaction ma¬ 
trix element in Eq. (38). The CM phase (if indeed 
exists) should be separated from the I phase by a 
CM-I quantum phase transition fixed point which 
is, however, located well beyond the limit of our 
one-loop RG. 


• As in other symmetry classes, we expect existence 
of a fixed point at 7 c = —00 and t ~ 1 (region 
marked by ‘SIT’) such that the transition between 
superconductor and insulator occurs through the 
separatrix connecting this fixed point and the triv¬ 
ial fixed point at t = 7 c = 0. At small values of t 
and 7 c < 0 the separatrix is parametrized by the 
following equation, t = 47 c [55] . 


• There is stable non-trivial fixed point at 7 * = 1/2 
and t* = 2/3 describing the critical metallic (CM) 


The RG flow (and the corresponding phase diagram) 
for equations (37) - (38) is shown in Fig. As in the 
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FIG. 9: (Color online) The case of broken spin rotational 
symmetry with Coulomb interaction, 7 s = —1: the RG flow 
obtained from numerical solution of Eqs. (371 and (38l. The 
arrows indicate the flow towards the infrared. The gray re¬ 
gion indicates the part of the phase diagram which is not 
accessible within one-loop RG equations. The lines describ¬ 
ing the flow to the superconducting (SC), insulating (I), and 
critical-metal (CM) phases are shown in red, green, and blue, 
correspondingly. Gray flow lines corresponds to the region of 
superconductor-insulator transition (SIT). 


case of preserved spin rotational symmetry, we stop the 
RG flow when either jyd reaches the value 1/t ^ I at 
a certain scale Lx (superconducting phase, red RG flow 
lines) or the resistance t reaches the value unity (insu¬ 
lating phase, green flow lines). In addition, we have now 
the critical metal phase (blue flow lines). 

The dependence of t on the length scale across the 


consecutive SG-TGM transitions is shown in Fig. 10 As 
one can see, there is a very narrow interval of to values 
in which the insulating phase separating the SG and GM 
phases exists. We mention that at not too large length 
scales (or, equivalently, at not too low temperatures) the 
resistance curves for SC, I, and CM phases cross each 
other. As expected, the Coulomb interaction suppresses 
the superconductivity, so that Lx > L^^^, and conse¬ 
quently, Tc < (as shown in the inset of Fig. [lol ). 

From Eq. (27) with n = 0 we find the following one- 


loop RG result for Zi^ in the case of Coulomb interaction: 


din Zi^ 
dy 




-b 27c -b 47) 


(39) 


As explained below Eq. (28), the fixed-point value of the 


C function d In Z^jdy determines the dynamical exponent 
controlling the temperature dependence of the specific 
heat. The one-loop result Eq. (39) yields for the critical 
metal Cy ~ -^vith z = 7/3. Since z > 2, the CM 
phase is stable with respect to the deviations of 7 ^ from 

Is = -1- 


FIG. 10: (Color online) The case of broken spin rotational 
symmetry with Coulomb interaction, 7 s = —1: dependence 
of t on the length scale across the transitions between su¬ 
perconducting (SC, red curves), insulating (I, green curves), 
and critical-metal (CM, blue curves) phases. The curves are 
obtained from numerical solution of Eqs. ( |37[ ) and ( |38[ ) for 
to = 0.03, 0.05, 0.058039, 0.0580414, 0.0580416, 0.0580418, 
0.0580419, 0.0580421, 0.07 (from bottom to top). With 
further increasing the Drude resistivity to, the system en¬ 
ters again the insulating phase (not shown here), see Fig. ^ 
The inset: dependence of on to- The value of the 

Cooper-channel attraction is 7co = —0.12. 


2. Short-ranged interaction 


In the case of short-ranged interaction RG equations 


(23) - (26) with n = 0 take the form 


^ -\+ fils) -Ic 


dt 
dy 

^ =-^(l + 7 s)( 7 s+ 27^ + 47^), 


(40) 

(41) 


dy ^'= + 2 1 


-(1 + 7c)7s + 27,'(1-27c) • (42) 


These RG equations are richer than Eqs. (33)-(36 1 and 
describe the RG flow in a three-dimensional space of t, 
Js, and 7 c: 


• There is a line of clean-Fermi-liquid fixed points at 
t = 7 c = 0 and arbitrarily 7 ^. The peculiarity of 
the present symmetry class is that the clean non¬ 
interacting fixed point t = 7 c = 7 s = 0 is attractive. 
It corresponds to a supermetal (SM) phase. 


• The line of fixed points at f = 0 and 7 c = —00 
corresponds to the superconducting phase. 


• As in all other symmetry classes, there is the insu¬ 
lating (I) phase with t = 00 not reachable within 
the one-loop RG equations. 
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FIG. 11: (Color online) The case of broken spin rotational 
symmetry with short-ranged interaction: a projection of the 
phase diagram on the plane 7co — 7 s 0 . The color indicates 
the value of to at which the QPT from SC to SM or I occurs. 
In regions 1, 2 and 3 the increase of to drives a sequence of 
QPTs: in 1 - SC-SM-I, in 2 - SC-SM-CM-I, and in 3 - SC-I- 
SM-I. The figure is obtained from numerical solutions of RG 
equations (401 - (421. 


• In addition to the critical metal phase at 7 * = 1/2, 
7 / = —1, and t* = 2/3 there is a second fully at¬ 
tractive fixed point at intermediate resistivity (i.e. 
on the border of applicability of one-loop RG equa¬ 
tions): 7 /* = —1/2, 7 ** = 0, and t** = 1. So, 
one-loop RG equations suggest a possibility of two 
different CM phases. 


• As in other symmetry classes, we expect existence 
of a fixed point t ~ 1 controlling the transition be¬ 
tween superconductor and insulator phases. Fur¬ 
ther intermediate-coupling (t ^ 1 ) fixed points 
control other emerging quantum phase transitions 
(SM-I, SC-SM, SM-CM, and CM-I). 


The phase diagram expected on the basis of RG equa¬ 
tions (40) - (42) is shown in Fig. 11 For 7 co < 0 the 
superconducting phase exists at small values of to ■ With 
increase of to the QPT to insulator or supermetal occurs 
for given values of 7 co and 7 so- With increase of 7^0 the 
superconducting phase proliferates. In the most part of 
the phase diagram the transition between superconduc¬ 
tor and insulator occurs at to 1- Interestingly, there 
is a region of the phase diagram in which a sequence 
of quantum phase transitions, SC-SM-I, SC-SM-CM-I or 
SC-I-SM-I, occurs as to grows (see Fig. 11). The dashed 
curve separating the region with multiple quantum phase 
transitions is parametrized by the condition 7^0 = 27 ^ 0 , 
see Appendix [B| 


A typical dependence of t on the length scale L across 
the SC-SM transition governed by increase of disorder 
(to) is illustrated in Fig. 12 for some initial values of 7so 
and 7co. As shown in the inset, in this case (dot-dashed 
line) the disorder suppresses the transition temperature. 
However, with increase of 7^0 one finds a nonmonotonous 
dependence of on to: at weak disorder is reduced in 
comparison with whereas at intermediate disorder 

Tc is larger than . There is the region in the phase 

diagram with small values of 7co < 0 and 7so < 0 in 
which T(, > (or more precisely. Lx < The 

enhancement of Tc in a certain range of bare couplings 
is in agreement with conclusions of our work m where 
renormalization group equations (40) - (42) with the right 
hand sides expanded to the lowest nontrivial order in 
7 s and 7 c were analyzed. In fact, for attraction in the 
particle-hole channel, 7^0 > 0 , a significant part of the 
phase diagram is occupied by the superconducting phase 
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with Tc > T^^^, see Fig. 

The mechanism of enhancement of the transition tem¬ 
perature is similar to one for the case of preserved spin 
rotational symmetry. For small initial values of 7 co and 
7 s 0 the renormalization of the Cooper interaction ampli¬ 
tude occurs in two distinct steps. At the first step of RG 
flow, the interaction is renormalized due to the presence 
of disorder (the terms proportional to t), while at the 
second step the standard BCS-type renormalization (the 
term 27 /) takes place. In the case | 7 so|, ITcoI «C to ^ I at 
the first step of renormalization we can linearize the RG 


equations (|40[)-(42 ) in interaction parameters and neglect 
the term — 27 /. Then, for 7 co < 7so/2 in the course of 
RG flow, the interaction amplitudes approach the BCS 
line 7 s = — 7 c, thus converting the repulsion in the singlet 
particle-hole channel into attraction. This is the conse¬ 
quence of the (weak) multifractality of the noninteracting 
fixed point. 

At some scale Li such that InLi/l = ^jt — 2!to 
the interaction couplings become of the order of the 
resistance: |7s,c| t- The resistance at this scale is 
t{Li) ~ [to(7sO — 27 co)]^^^ ^ to- We note that if the 
length scale Lx is reached before Li which is typical 
foi' to ^ niax{|7co|, |7so|} one can find suppression of 
the transition temperature instead of enhancement (see 
solid curve in the inset to Fig. 12). After the length 
scale Li the second step of RG flow starts, where in 
Eq. ( [ 4 ^ one can neglect terms proportional to t com¬ 
pared to the disorder-independent term — 27 /. Thus, the 
Cooper interaction 7 ^ flows according to the standard 
BCS RG equation for a clean system with the initial 
value 7 c(Ti) ^ t{Li) rather than 'jco- Hence, we find the 
following rough estimate for the transition temperature: 
lnl/(TcT) -- 21nLi/l-(-l/|7c(^)| -- [to{lsO-‘2lco)V^‘^ > 
lnl/(T/^C"5r), see Appendix for details. 

In full analogy with the case of preserved spin-rotation 
invariance, if the disorder scattering rate l/r exceeds 
the Debye frequency wu, the starting point of the RG 
flow will be likely located not far from the BCS line, 
7 s 0 = — 7*0 = — 7 c 0 (see Sec. HC). For such initial condi- 
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FIG. 12: (Color online) The case of broken spin rotational 
symmetry with short-ranged interaction: dependence of t on 
the length scale across the QPT between superconducting and 
supermetallic phases. The curves are obtained from numerical 
solutions of RG equations ( [40| ) - (421 for jco = —0.04, 7^0 = 
—0.1, and to = 0.05, 0.07, 0.1, 0.12, 0.15, 0.2 (from bottom 
to top). Inset: Dependences of Tc/{Tc + on to for 

7co = —0.04 and 7^0 = —0.005 (solid curve), for 7co = — 7 s 0 = 
—0.05 (dashed curve) and for 7co = —0.04 and 7 so = —0.1 
(dot-dashed curve) are shown in the inset. 


tions, the dependence of Tc/(Tc -I- on to is shown 

in the inset to Fig. [I^ by dashed curve. As in the case 
of preserved spin invariance, for bare interaction param¬ 
eters on the BCS line, there is no initial decrease of tran¬ 
sition temperature with increase of to- For a given jco, a 
deviation from the BCS line in the initial conditions to¬ 
wards larger (smaller) values of 7so increases (decreases) 
the relative enhancement of the transition temperature, 
Tc/T^^^. We emphasize that inspite of the antilocaliza¬ 
tion at the noninteracting fixed point the multifractality 
enhances the interaction in the Cooper channel. 


V. RESISTANCE IN ZERO MAGNETIC FIELD 


Within the NLSM approach, physical observables can 
be written as correlation functions of the matrix field Q. 
In particular, the conductivity obtained by evaluating a 
linear response to an electromagnetic field in the frame¬ 
work of the NLSM theory with the action ([^ can be 
expressed in the following way: 

a{iuJn) =- Y|^^Tr[J“,Q(r)][J“„,Q(r)]^ 

+ ^J dr-'(Tr J“Q(r)VQ(r) 

xTrJ“„Q(r')VQ(r')). (43) 



FIG. 13: (Color online) The case of broken spin rotational 
symmetry with short-ranged interaction: the color plot for 
the ratio of Tc/{Tc -|- for to = 0.11. The dashed curve 

separates the regions with Tc < and with Tc > . 

The solid line indicate the boundary between superconduc¬ 
tor and supermetal. The figure is obtained from numerical 
solutions of RG Eqs. (401 - (42l. 


Here = 2TTTn is a Matsubara frequency, expectation 
values are defined with respect to the action 0 , and 


ja _ ^30 ^00 JO 


^30 + too 


(44) 


As usual, the static conductivity cr can be obtained after 
the analytic continuation of Eq. (431 to real frequencies: 
iojn —> w -|- iO^ and, then, taking the limit a; —>■ 0. At the 
classical level, Q — A, one finds a = g. 

The RG equations derived in this work describe renor¬ 
malization of the couplings in the NLSM action with the 
running spatial scale L. As such, these equations yield 
physical observables (e.g., the resistivity) of a finite size 
sample at T = 0. At finite temperature T, conductivity 
can be evaluated in two steps. At the first step, the ac¬ 
tion ([^ is renormalized from the energy scale 1/r down 
to T. Thus the bare parameters in the action 0 are 
substituted by the parameters at the length scale Lt'- 
g 2/[7rt(Lr)], ls,t,c ls,t,c{LT), and -J- Z^{Lt). 
Their dependence on Lt is governed by RG equations 
(27). At the second step, the Kubo formula (43) 


is evaluated under assumption that in the NLSM action 
([^ Q fields are restricted by the temperature in the ul¬ 
traviolet. We assume Lt to be such that t[LT) <C 1 and 
|7 c(Lt)| ^ 1- Then, the conductivity can be written as 


cr(T) ~ ^-7 c(L 7') In — 

^ ' TTt{LT) 2 Lt 


(45) 


This result illustrates the fact that at finite tempera¬ 
ture there is always a difference between the physical 
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FIG. 14: (Color online) Dependence of the resistance p on temperature T and Tc on to for the cases of Coulomb [(a), (b) and 
(c)] and short-ranged [(d), (e) and (f)] interaction. Parameters are the same as in Figs. [6 10 ^ and [l^ respectively. 
Figures (a), (b), (d) and (e) correspond to the case of preserved spin rotational symmetry, (^ and (f) - broken spin rotational 
symmetry calculated from Eq. | |45[ ). The dashed parts of the curves correspond to the “fluctuation” region T < Tx- Dotted 
red curves in panels (a) and (d) show temperature dependences of t for T > Tx- The dependence of Tc on to is shown in the 
insets (solid curve). The dependence of Tx on to is shown by dashed line in the inset to figure (a). For all other insets Tx 
coincides with Tc within our accuracy. 


resistance p = 1/a and the coupling parameter t in the 
NLSM action. The second term in Eq. ( |4^ is the Maki- 
Thompson contribution which is the dominant 

inelastic contribution to the resistivity. (We neglect the 
smaller Aslamazov-Larkin contribution HD].) 

Comparing two terms in Eq. (45), we find that the 


inelastic contribution becomes of the same order as the 
renormalized impurity-scattering conductivity (given by 
the NLSM coupling) at the scale Lx determined by the 
condition t{L)\'yc{L)\ = 1. Remarkably, this is the same 
scale as is found from the condition of validity of the 
one-loop RG, see Sec. |III[ Thus, at scales shorter than 
Lx, or, equivalently, at temperatures larger than Tx, 
the temperature dependence of the physical resistivity is 
dominated by the RG behavior of t{L). Upon approach¬ 
ing the transition temperature, the role of the inelastic 
contribution controlled by 7 c(A) increases. In the nar¬ 
row temperature interval < T < Tx the conductivity 
is dominated by the inelastic contribution. As we have 
already discussed in Sec. |IV[ the width of this interval is 

(46) 


{Tx-Tc)/Tc^t{Tx). 


0 . 

In fact, the behaviour of the conductivity in the fluctu¬ 
ation region, \T — Tc\/Tc ~ t{Tx), is additionally affected 
by superconducting phase fluctuations. These fluctua¬ 
tions lead to the BKT character of the actual supercon¬ 
ducting phase transition. The corresponding shift of the 
transition temperature is, however, small (as has been 
already mentioned in the end of Sec. IV A 11 and is not 


At the transition point the running coupling di¬ 
verges, | 7 c(Tc)| = oo and Eq. (451 formally yields p(Tc) = 


important for our results. More specifically, as was ar¬ 
gued by Beasley, Mooij, and Orlando [S^, the shift is 
determined by disorder strength, (Tc-T-q-xt)/T^ ~ t, for 
t 1, see also Refs. |84l85j . The behavior of the true 
p(T) in the fluctuation (BKT) region will be addressed 
elsewhere [54|. Let us only mention here that the value 
of the NLSM coupling t{Tx) at the entrance to the fluc¬ 
tuating region determines [M] the stiffness of the phase 
fluctuations in the BKT region Tbkt <T <Tx and the 
relevant value of t in the shift of Tbkt with respect to 
Tc- Because of the renormalization of t, this value may 
strongly differ from the bare (high temperature) Drude 
value of the resistance to (cf. Ref. |14j). 

Ignoring the above-mentioned subtleties of the behav¬ 
ior of p{T) in the narrow fluctuation region |T —Td/T^ ^ 
t{Tx) around the mean-field Tc, the behavior of the elec- 
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trical resistance in the whole range of temperatures is well 
described by p{T) = l/a{T) with a{T) given by Eq. (45). 
To relate the length scale with the transition temper¬ 
ature Tc one needs to know the relevant value Dc of the 
diffusive coefficient D ~ It is obtained by using 

the values for t and at the scale Lx (which is the bor¬ 
der of validity of the one-loop RG and simultaneously is 
the beginning of the fluctuation region). The transition 
temperature is given by Tc ~ DcL~‘^, see Eq. (291. 

We present the dependence of p{T) obtained in ac¬ 
cordance with Eq. (451 (for simplicity, we dropped the 
logarithmic factor in the inelastic term; this does not 
qualitatively affect the plots) for the cases of preserved 
and broken spin rotational symmetry and Coulomb and 
short-ranged interactions in Fig. We mention that 
due to Maki-Thompson correction the resistance drops 
very fast (dashed curves in Fig. 14) since \^c{Lx)\ ^ 1 
and, consequently, {Tc — Tx)/Tx ~ p{Tx) ^ 1- There¬ 
fore, the temperature Tx, at which resistivity has the 
maximum, can be used as an estimate of the supercon¬ 
ducting transition temperature. In order to plot p{T) in 
the fluctuation region T < Tx where the one-loop RG 
becomes insufficient, we evaluated ^c{T) near Tc keeping 
only the BCS term — 27 ^ in the RG equation for 7 ^. 


We note that our approach for evaluating temperature 
dependence of resistivity is different form that of Refs. 
[TSTTS] . In these works the full set of the first-order per¬ 
turbative quantum corrections to the conductivity was 
computed at finite temperature in the presence magnetic 
field. However, such approach assumes that T and H de¬ 
pendent corrections to the bare Drude conductivity 1/to 
are small. In our approach we split effects leading to 
temperature variation of the conductivity into two parts: 
those related to virtual and real processes. Then virtual 
processes are taken into account within the RG formal¬ 
ism (they are included in the renormalisation of t). This 
allows us to consider also situations with strong renor¬ 
malisation, t{LT) S> toj i- 6 -; cases in which “quantum 
corrections” are large in comparison with 1 /to- 


VI. MAGNETORESISTANCE 

A transverse magnetic field introduces an additional 
length scale Ih (magnetic length) into the problem. In 
what follows we assume that it is larger than the mean 
free path, I < Ih- Let us start from T = 0. In case 
of weak magnetic field, Ih ^ Ac, the superconducting 
instability at L = Lc remains unaffected. For strong 
magnetic fields, Ih ^ Lc, the superconducting phase is 
destroyed by magnetic field since the growth of | 7 c| within 
RG equations is stopped at the length scale Ih- Then the 
critical magnetic field can be estimated as Ih^ = Lc- This 
results in the standard relation, He ~ Tc/Dc- At H < 
He the physical resistance p should vanish in the infinite 
system, L = 00 - For systems of finite size L > Lc, the 
resistance is not zero and is determined by the nontrivial 
configurations of the order parameter in the presence of 


magnetic field. 

At finite temperature, the critical magnetic field is a 
function of temperature. However, at T <C T/ this ef¬ 
fect is small and can be neglected. In order to evaluate 
the magnetoresistance at T ^ Tc and for H > He we 
use a two-step RG procedure. At the first step of RG 
the interaction in the Cooper channel grows towards in¬ 
stability. The first step ends at the length scale Ih with 
some values t{lH) and 'js,t,c{lH)- At L > Ih the cooperon 
modes become ineffective and the second step of the RG 
procedure starts. For length scales Ih < L < Lt RG 
equations do not contain cooperon contributions: 


dt 
dy 
djs 
dy 
djt 
dy 
din Z,^ 
dy 


fils) +nf{jt) 


“ 2(1 + 7s)(7s +«7t)> 
-^(1 + 7t) Is - in-2)-ft 


(7s + - 


(47) 

(48) 

(49) 

(50) 


Here y = \tiL/Ih and the initial values of couplings in 
Eqs. (47) - ( [^ are given by tiln) and 7 s.t,c(^ff)- Sim¬ 
ilarly to the case of zero magnetic field, in order to find 
the physical resistance p{T, H) one needs to take into ac¬ 
count non-RG corrections to conductivity due to super¬ 
conducting fluctuations. In particular, for | 7 c(^//)| >> 1 
the leading-oder correction (which is the dominant one 
for magnetic fields not too close to the critical field He) 
[77] yields: 


a{T,H) = 


(51) 


It should be emphasized that tiLx) in Eq. (51) in fact 


depends on the magnetic field via the two-step RG pro¬ 
cedure. At H = He the physical resistivity p{T, H) 
should vanish. This happens due to higher-order Cooper- 
channel corrections of non-RG type in Eq. (51) which be¬ 


come important in the fluctuation region (i.e., near He) 
and make p{T,Hc) = 0 in spite of finite value t{Lj’)- 
Let us mention that, due to renormalization on bal¬ 
listic scales, L < I, |78|, one can also expect an effect 
of magnetic field on the initial values of parameters (to, 
7 c 0 : Im and 7 so) for RG equations (23) - ^7\- This leads 


to a shift of the length scale at which 7 ^ diverges. As a 
consequence, an additional dependence of transition tem¬ 
perature on H appears. We do not take this effect into 
account. 

The dependence of the resistivity on perpendicular 
magnetic field H at different temperatures below Tc is 
illustrated in Fig. The H dependence of resistiv¬ 
ity at fixed T < Tc is qualitatively similar to the p(T) 
dependence at zero magnetic field. For all four symme¬ 
try classes the magnetoresistance shows a maximum that 
grows with decreasing temperature. It is worth mention¬ 
ing that the maximum should become arbitrarily high 
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FIG. 15: Dependence of the resistance p on perpendicular 
magnetic field H for the cases of preserved [(a) and (b)] and 
broken [(c) and (d)] spin rotational symmetries. Figures (a) 
and (c) correspond to the case of Coulomb interaction, (c) 
and (d) - short-ranged interaction. Solid curves are obtained 
in from Eq. (511. The dashed lines indicate the drop of 
magnetoresistance to zero at H = He- The parameters used 
are as follows (a) 7 s = —1, 7 c 0 = —0.45, 7 t = 1, to = 0.1, 
T = Tc/2, Tc/4, Tc/ 8 , (b) 7 so = 0.1, 7 co = —0.1, 7 to = —0.1, 
to = 0.05, T = Tc/2, Ts/4, n/16 (c) 7 s = -1, 7 c 0 = -0.12, 
to = 0.053, T = Ts/2, Tc/4, T,/8, and (d) 7 so = 0.05, 7^0 = 
-0.05, to = 0.2, T = n/2, n/16, rc/256. 


as T —0, yielding a giant magnetoresistance in agree¬ 
ment with experimental observations. We do not plot 
these curves here since we cannot controllably evaluate 
the resistance for t >1 within the one-loop RG. 

In addition to Ih, the perpendicular magnetic field in¬ 
duces another length scale Iz related to the Zeeman split¬ 
ting. Usually one expects that Ih Iz since the lat- 
ter can be estimated as tz /\/(l + jt{lz))t{lz)gL- 
Here gz stands for the Lande g-factor. Due to Zeeman 
splitting the magnetic field suppresses the triplet diffu- 
son modes with = ±1 and singlet and triplet cooperon 
modes with S'z = 0 at length scales L > Iz- 

In the case of fully broken spin-rotational symmetry 
there is no triplet modes, and thus the orbital and Zee- 
man effects of magnetic field on RG equations are the 
same. Therefore, there is in fact only one length scale 
associated with magnetic field Ihz = inin{Z//, ^ 2 }. A 
similar conclusion holds for the case of partially broken 
spin-rotational symmetry, n = 1. Therefore, in the ab¬ 
sence of spin-rotational symmetry we do not need to con¬ 
sider the effect of Zeeman splitting separately. 


For the case of preserved spin-rotational symmetry, the 
orbital and Zeeman effects are different. For Ih ^ Iz, the 
RG equations (47) - (50) should be modified at length 
scales L > Iz since two out of three triplet diffusive 
modes becomes massive and do not lead to infrared di- 




FIG. 16: Dependence of the resistance p on parallel magnetic 
field H for the case of preserved spin rotational symmetry 
for the case of Coulomb (a) and short-ranged (b) interaction. 
Solid curves are obtained from Eq. (56 1 . The dashed lines 
indicate the drop of magnetoresistance to zero at H = Hez- 
The parameters used are as follows (a) 73 = —1, 7 co = —0.45, 
7 t = 1, fo = 0.1, T = Tc/2, Tc/4, T^/S, (b) 7so = 0.1, 
7 co = -0.1, 7 to = -0.1, to = 0.05, T = n/2, n/A, T 3 /I 6 


vergences: 


dt 

dy 

djs 

dy 

djt 

dy 

din 

dy 


f{’ls) +fht) , 


t 

2 

t 

2 


(1 + 7s)(7s +7t): 


(1 + 7t) 


'Is + 7t 



(52) 

(53) 

(54) 

(55) 


Here y = InL/lz and the initial values of couplings in 
Eqs. (52) - (55) are given by t{lz) and 7 s,t,c(^z)- There¬ 
fore, in this case, a three-step RG scenario is realized. 

In the opposite case, Ih ^ Iz, the system at length 
scales I < L < Iz is described by RG equations (23) - 


(27) (with n = 3). Then for Iz < L < Ih RG equations 


transforms into Eqs. (521 - ([55| with /( 7 s) substituted by 


1 -h /{"fs) (weak-localization correction remains intact in 
the presence of Zeeman splitting only). For larger length 
scales, L > Ih, the system is governed by RG Eqs. (52) 


(55). The suppression of weak-localization correction 


does not change qualitative behavior of the RG flow. 

In the case of a parallel magnetic field, the scale Ih 
does not appear since the orbital effect of magnetic field 
can be neglected provided Iz ^ Lt ^ ^i/d, where d is 
the typical width of the film. We thus obtain a two-step 
RG scenario in which RG equations are modified at the 
length scale Iz- As before, in order to find the physical 
resistance p(T, H) one needs to take into account non- 
RG corrections to conductivity due to superconducting 
fluctuations. For |7c(^ff)| 1 we get [75] 


aiT,H) = 


TtiLr) 


ln|7c(^z)|- 


(56) 


The parallel-field magnetoresistance for the case of pre¬ 
served spin-rotational symmetry is illustrated in Fig. 
for several values of temperature below T^. The depen¬ 
dence of resistivity on parallel magnetic field at fixed 
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temperate below Tc is essentially different from the p{H) 
dependence in the case of transverse field. First, the max¬ 
imum at an intermediate field is much less pronounced in 
the case of parallel field. Second, the parallel-field resis¬ 
tivity increases with H in strong fields, contrary to the 
case of transverse magnetic field. 


VII. DISCUSSION 

In this section, we discuss our results and their impli¬ 
cations (in particular, for the phase diagrams of SITs), 
relation to previous works, as well as limitations and pos¬ 
sible generalizations of the RG scheme used. 


A. Relevant superconducting systems 

1. Symmetry of the order parameter 

First of all, we note that our theory is derived for con¬ 
ventional (BCS) s-wave superconductors, where the ef¬ 
fect of s-wave non-magnetic impurities on the supercon¬ 
ducting gap and is absent at the semiclassical level 
(“Anderson theorem”). Our theory can be generalized 
to describe the multiband case with s-wave (or s±-wave) 
pairing. On the other hand, in unconventional (p-wave 
or d-wave) superconductors, impurities do suppress the 
superconductivity. As a result, the diffusion regime does 
not develop there: either the pairing is so strong that the 
superconductivity is established already on the ballistic 
scales, or disorder kills the superconductivity. Therefore, 
in such systems the enhancement of superconductivity 
by localization (which occurs on the diffusive scales in 
s-wave superconductors) is impossible. 

However, in such superconductors a secondary super¬ 
conducting transition due to the pairing of Dirac quasi¬ 
particles is possible (which may change the true gap sym¬ 
metry as, e.g., d —)■ d+is, thus opening the superconduct¬ 
ing gap at the nodal points of the spectrum), see e.g. Ref. 
m for review. This transition can be described by the 
RG equations generalized for the novel symmetry classes 
(see Refs. mm)- Furthermore, the peculiar form of the 
Fermi surface near half filling (nesting) may lead to addi¬ 
tional emergent symmetries specific to this problem [52] . 
In particular, various novel interaction couplings would 
be possible by the enhanced symmetry. Importantly, the 
Coulomb interaction between the quasiparticles in this 
system is screened by the d-wave condensate. Thus, one 
can expect a disorder-induced enhancement of the criti¬ 
cal temperature for the secondary superconducting tran¬ 
sition. 


2. Macroscopic homogeneity vs granularity 

In this paper we assume that the system is macro- 
scopically homogeneous and do not discuss granulated 


superconductors characterized by weak (Josephson) tun¬ 
neling between macroscopic superconducting islands. In 
granular systems, additional energy scales appear such as 
Josephson and charging energies. We expect, however, 
that the peculiarities of inhomogeneous superconduc¬ 
tors, while leading to emergence of intermediate crossover 
regimes, do not affect the universality of the (zero-T) SIT 
governed by the symmetries of the system. The situation 
resembles the problem of Anderson metal-insulator tran¬ 
sition which is believed to be universal independently 
of whether the microscopic disorder model is “homo¬ 
geneous” (e.g. white-noise disorder) or inhomogeneous 
(tunnel-coupled grains). 

At the same time, the finite-T behavior of the resistiv¬ 
ity in granular systems will be influenced by the presence 
of additional energy scales and thus differ from that of a 
homogeneous system. On the other hand, the behaviour 
close to the transition will be governed by a similar BKT 
physics both for granular and homogeneous systems, see 
the discussion in SecjVj 


B. Screening of long-ranged Coulomb interaction 


Above, we have considered separately the two models 
of electron-electron interaction: long-range Coulomb in¬ 
teraction and short-ranged interaction. In the latter case, 
the superconductivity was shown to be enhanced by An¬ 
derson localization in a wide parametric range. In realis¬ 
tic electronic systems, there are two mechanisms that can 
suppress the Coulomb interaction and make it effectively 
short-range in a certain interval of length scales: (i) large 
dielectric constant of the medium, and (ii) screening by a 
nearby external metallic layer which results in a less sin¬ 
gular dipole-dipole type interaction at scales larger than 
the distance to the gate. 

In the presence of a dielectric medium, the interaction 
constant in the singlet channel acquires a momentum de¬ 
pendence: 


is{q) = Is 


(1 + 7s) 


>t 

+ g’ 


27re^ dn 
e dp 


(57) 


Here 7 s is the irreducible short-ranged part of the singlet 
interaction amplitude, e is the dielectric constant of the 
medium, >f is the inverse screening length, and dn/dp is 
the thermodynamic density of states (which is not renor¬ 
malized by the interplay of disorder and interaction). 

Usually, the condition >c~^ ^ I (the screening radius 
is smaller than the mean free path) is fulfilled and for 
length scales L ^ I one finds 7 s = —1 which is a hallmark 
of long-ranged Coulomb interaction. However, for large 
dielectric constant the opposite relation, I ^ is pos¬ 
sible. In this case, at length scales I ^ L ^ the long- 
ranged Coulomb interaction provides small contribution 
to 7 s and is indistinguishable from the short-ranged inter¬ 
action within the RG. If the scale Lc is smaller than >c~^, 
then the long-ranged Coulomb interaction does not affect 
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the transition temperature. Therefore, for large dielec¬ 
tric constant such that <C 1, the long-range Coulomb 
interaction does not influence the superconducting tem¬ 
perature provided the following condition holds: 


Tc. 


(58) 


When condition (581 is not fulfilled, the long-ranged 


nature of Coulomb interaction screened by high dielectric 
constant becomes effective at large scales L > before 
the superconductivity occurs. While for scales shorter 
than the screening radius the coupling constant 7 c is en¬ 
hanced by short-range interaction as compared to the 
BCS result, at larger scales the Coulomb repulsion starts 
working in the opposite direction. As a result, one en¬ 
counters the competition between the enhancement and 
suppression of the superconductivity. In this situation, a 
more general scheme of including Coulomb repulsion is 
necessary. 

The simplest generalization of the RG procedure would 
be then a two-step RG. At the first step, for L < 
one uses the short-ranged RG with the initial values of 
all interaction couplings determined by the short-range 
attraction (BCS line). At the second step, for L > 
the RG equations are switched to the Goulomb case with 
7 s = —1 and the initial values of other couplings given 
by the outcome of the first step. However, within this 
two-step procedure the singlet amplitude 7 s is instantly 
switched at L ^ from the value dominated by the 
phonon-induced attraction, 7 s(A), to the Coulomb dom¬ 
inated value 7 s = —1, implying the change of its sign. 

In order to smoothly describe the crossover regime, an 
interpolating flow equation for the coupling 7 s defined in 
Eq. (571 can be derived by replacing the momentum by 
. In particular, for the case of preserved time and 


spin-rotational symmetries this yields [cf. Eq. (191] 


^ = -2(l+A^)(7^ + 37*+27c + 47e^)^^ 


- (l+7s) 




+ 1 


(59) 


In this flow equation we have introduced the new coupling 
Z]^ = kL satisfying 


dZi^ 

dy 


= Zl. 


(60) 


In Eq. (601 we have used the fact that k is not renor¬ 


malized by interactions, since it is determined by elec¬ 
tron charge and the thermodynamic density of states 
dn/dy, [52]. If the background medium is characterized 
by a momentum-dependent dielectric function e{q), this 
would modify Eq. (60) accordingly. The RG flow gov¬ 


erned by Eqs. (59) and (60) can be viewed as a two- 


step RG procedure with a short-ranged singlet amplitude 
7 s ~ 7 s(l -I- Zl) — Zl at the first step. 

We plot the results of the numerical evaluation of the 
renormalization of the Cooper-channel coupling 7 c for 



FIG. 17: (Color online) Effect of a high dielectric constant 
on the renormalization of Cooper-channel coupling ■jc- Black 
curve corresponds to solution of RG Eqs. (33) - (361 for 7co = 
7*0 = —730 = —0.1, and to = 0.05. Blue curve corresponds 
to the behavior of 7 c on L for clean BCS case. Red curve 
is obtained from the numerical solution of the two-step RG 
equations for ln(l/xrZ) = 2. Red dashed curve is obtained 
from the numerical solution of the crossover RG Eqs. | |59[ ) 
and (60l. Inset: Dependence of Tc/{Tc + on ln(l/xrl) 

for two-step (solid) and crossover (dashed) RG equations. 


ln(l/>fZ) = 2 using the two-step and interpolating RG 
procedures in Fig.fTTj One can see that the enhancement 
of the superconductivity at the first (short-ranged) step 
of the RG is more important than the suppression at the 
second (long-ranged) step. As shown in the inset, the 
overall enhancement of the superconductivity takes place 
(for chosen values of the bare interactions and resistivity) 
for > 3 -b 4. For larger screening lengths both 

the two-step and interpolating RG procedures yield close 
results for the enhancement of Tc- 

The long-range Coulomb repulsion can also be screened 
by a nearby metallic layer. Specifically, the electron- 
electron repulsion can be considered as short-ranged on 
scales L larger than the spacer width Wg- When the 
mean-free path is larger than Wg-, we have the short- 
range case from the very beginning. In the opposite case 
Wg ^ I, without additional screening by the dielectric 
medium (i.e. for kL 3> 1) the RG procedure corre¬ 
sponds to the Goulomb case 7 s = —1 (up to small cor¬ 
rections of the order of {KWg)~^ <C 1) for scales L < Wg. 
For larger scales the interaction becomes of the dipole- 
dipole type, but the singlet interaction constant inher¬ 
ited from the first step remains Coulomb-like, 7 s — — 1. 
Therefore, a metallic layer placed at the distance Wg » I 
is not sufficient to screen the Coulomb repulsion such 
that the superconductivity would be enhanced. How¬ 
ever, the combination of the screening by a medium with 
large dielectric constant (see above) and by the metal¬ 
lic layer does lead to the enhancement of the supercon- 
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ductivity as compared to the cases when these screening 
mechanisms are considered separately. Indeed, these two 
mechanisms make the interaction effectively short-ranged 
(with | 7 s| < 1) for short and large scales, respectively. In 
particular, for xws ^ 1 there is no room for the Coulomb 
regime at all. 


C. Enhancement of superconductivity for 
short-range repulsion 

Most of experiments on the superconducting transition 
in 2D films have been performed without screening the 
long-range component of the interaction. It is desirable 
to explore whether the mechanism of the enhancement of 
superconductivity addressed in the present work may be 
employed in practice to obtain structures with substan¬ 
tially enhanced T^. The key condition is a suppression 
of the long-range component of the Coulomb interaction 
[S7155] . This opens a new way for searching novel mate¬ 
rials exhibiting high-temperature superconductivity: one 
needs the combination of a large dielectric background 
constant and disorder in layered structures. 

As mentioned in Introduction, 2D superconductivity 
has been recently realized in interfaces between two ox¬ 
ides, in particular, in LaAlOa/SrTiOa interfaces [I6II7j . 
These systems possess unique electrostatic properties ow¬ 
ing to the giant dielectric constant of SrTiOa. In partic¬ 
ular, the long-range component of the Coulomb inter¬ 
action is expected to be strongly screened in such ma¬ 
terials. Although currently, the highest Tc reached in 
such materials is rather low as compared to high-T^ ma¬ 
terials, the dependence of on the conductivity of a 
normal state is non-monotonic, which agrees with the 
localization-induced mechanism of the superconductiv¬ 
ity enhancement. Further investigations are required to 
identify the ways for increasing in strongly screened 
oxide interfaces, and to analyze optimal materials, struc¬ 
ture design, and operation regimes, depending on the 
microscopic details. 

A possible route for increasing the superconducting 
transition temperature in these materials is based on fur¬ 
ther suppression of the long-range Coulomb interaction 
by designing a double-interface structure with a LaAlOa 
layer sandwiched between two SrTiOa oxides. In such 
a setup, already ten atomic layers of LaAlOa are suf¬ 
ficient, so that the screening properties of the sample 
would be fully determined by the giant dielectric constant 
of SrTiOa. At the same time, the two interfaces would 
be coupled by the interlayer interaction, similarly to the 
Coulomb drag problem in double-layer structures. The 
corresponding generalization of the sigma-model would 
include an additional degree of freedom (a pseudospin in 
the interface space). Furthermore, the doping of SrTiOa 
layers away from the interfaces can produce an effective 
metallic gate made of the same material. 

A simpler setup would involve an amorphous supercon¬ 
ducting film placed on a SrTiOa substrate (again possi- 



FIG. 18: (Color online) The case of broken spin rotational 
symmetry. Comparison of superconducting transition tem¬ 
perature for the Coulomb interaction {T^) and for the short- 
ranged interaction on the BCS line The color indicates 

the ratio {T^ + T^^^)/+ The dashed curve cor¬ 

responds to the boundary of SC phase in the case of Coulomb 
interaction. 


bly doped away from the interface) with high dielectric 
constant instead of more conventional Si 02 or AI 2 O 3 in¬ 
sulating substrates typically used in experiments on the 
SIT. On the other side of the substrate one can place a 
metallic gate, thus realizing both mechanisms of screen¬ 
ing discussed above. An interesting possibility of ar¬ 
ranging a closely located metallic layer is provided by 
BN-Graphene heterostructures [90] with gated graphene 
layer serving as a metallic gate and BN playing a role of a 
thin spacer. In this situation, a generalization of the RG 
equations (23) - (27) to the case of two layers (similar to 
Ref. [59]) needs to be done. 

Recently, the superconductivity has been studied in 
layered material LijjZrNGl [21]. It was found that with 
increase of doping level x the transition from insula¬ 
tor to superconductor occurs at a: « 0.05. The criti¬ 
cal resistance is close to /i/2e^. The temperature de¬ 
pendence of the resistivity, p{T), measured across the 
transition is qualitatively similar to one shown in Fig. 
[T4| 1. Near the SIT the superconducting transition tem¬ 
perature increases with decrease of doping level: from 
Tc{x « 0.12) « II K to Tc{x « 0.05) « 16 K. Such be¬ 
havior of superconducting temperature is in suit with de¬ 
pendence Tc{to) on BCS line predicted by our theory (see 
the inset to Fig. |^and Fig. [^. Similar nonmonotonous 
dependence of Tc on disorder was measured in W and Mo 
based films (for an overview, see Ref. [H]). 

In the case of repulsion in the particle-hole channel, the 
enhancement of superconductivity by Anderson localiza¬ 
tion (in comparison with the corresponding clean system) 
occurs in a certain range of (not too strong) interaction 
and (not too weak) disorder, see Figs. [^ [^ and [^ It 
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should be stressed, however, that the critical tempera¬ 
ture for short-range interaction is predicted to be always 
higher than for unscreened Coulomb interaction when 
other parameters are kept fixed, see e.g. the inset to 
Fig. 17 Therefore, we propose to perform benchmarking 
experiments, measuring Tc in the same superconducting 
film placed on the substrate with high dielectric constant 
(say, STO-material that screens long-range Coulomb re¬ 
pulsion) and on the reference substrate with a not too 
high e (say, Si 02 or AI 2 O 3 ). The experiments should 
be performed for sufficiently dirty samples (but still on 
superconducting side of the SIT), since the stronger dis¬ 
order leads to a stronger difference between the critical 
temperatures in the two cases, see Fig. [Tsl 


D. Relation to numerical results 

Recent numerical calculations [57155] demonstrate that 
disorder may indeed enhance the superconductivity in a 
certain range of parameters. These results should be con¬ 
trasted with numerical simulations of a two-dimensional 
disordered Hubbard model with strong on-site attraction 
in small-size systems that yielded a monotonous suppres¬ 
sion of Tc with increasing disorder [53]. The physics be¬ 
hind the results of our work develops for not too strong 
disorder and interaction, whereas Ref. [S3] focussed on 
the opposite limit. Specifically, we predict the enhance¬ 
ment of Tc by Anderson localization in 2D when both dis¬ 
order and interaction are weak: | 7 co| to ^ ■\/bA)T 
1. In terms of the disordered Hubbard model used in nu¬ 
merical simulations [89] this regime corresponds to the 
following range of parameters: |C/| ^ H ^ \/W\ ^ 1 
where U and V stand for dimensionless interaction and 
disorder. 

It is the strong interaction U which allowed the au¬ 
thors of Ref. [53] to extract the information on supercon¬ 
ducting properties from the simulation on a rather small 
system of 8 x 8 sites. As seen from Fig. [^ in the strong¬ 
coupling regime, our theory agrees with the numerics of 
Ref. [53]. Indeed, for strong attraction and strong disor¬ 
der our theory predicts a suppression of the mean-field 
Tc (and hence a suppression of the true critical tempera¬ 
ture). The point is that, in realistic systems, the attrac¬ 
tion is normally considerably weaker (otherwise Tc would 
be given by the Debye energy and no challenge of obtain¬ 
ing high-temperature superconductivity existed) and this 
is precisely the range of 7 c where the enhancement of su¬ 
perconductivity is expected according to our predictions. 
The results of Ref. [53] may in addition reflect the dif¬ 
ference between Tc and Tbkt in the strongly disordered 
case, see Sec. [^ If one extrapolates the Beasley-Mooij- 
Orlando estimate [83] to the regime of strong disorder, 
~ I, one gets {Tc — Tbkt) /7c I- This correlates 
with the numerical findings of Ref. [S3] : for strong disor¬ 
der Tbkt niay be significantly lower than the mean-field 
Tc (this difference might be important near the SIT). 

In fact, disordered Hubbard model contains all ingre¬ 


dients required for the enhancement of superconductivity 
by multifractality, but the range of optimal parameters 
requires large system sizes. To verify our prediction nu¬ 
merically within the disordered Hubbard model, one has 
to use weaker interaction and disorder and hence larger 
system sizes oi N x N sites, where due to the logarith¬ 
mic renormalization of couplings in 2D N depends expo¬ 
nentially on the inverse disorder strength. Rough esti¬ 
mates yield at least A^ 30 q- 50 for the minimal system 
size where the enhancement can be detected. Indeed, in 
Ref. |87| . where an enhancement of the superconductiv¬ 
ity by disorder in a honeycomb lattice was detected for 
a certain range of parameters, the number of sites in the 
attractive Hubbard model was 900^-1600, in consistency 
with the above estimate. At the same time, such system 
sizes are still much smaller than the sizes of real macro¬ 
scopic systems where the regime required for a strong 
enhancement of Tc by our mechanism can be realized. 


E. Beyond one-loop RG: Structure of the phase 
diagram 


In Sec. |IV| we have analyzed the one-loop RG equa¬ 
tions for the cases of preserved and broken spin-rotational 
symmetry. As we have shown, the one-loop precision 
is applicable for tmax{I, lyd} < I, and therefore some 
fixed points of the full phase diagram remain unaccessi- 
ble at this level, see gray areas in Figs, [^and]^ Here 
we discuss an expected structure of the full phase dia¬ 
gram, going beyond the one-loop RG equations. We will 
focus on the SIT part of the phase diagram, first disre¬ 
garding the complications related to the appearance of 
additional phases such as ferromagnetic (FM) phase for 
the case of preserved spin-rotational symmetry and the 
critical metal (GM) for the spin-orbit case. Furthermore, 
for simplicity, we consider the Goulomb case, 7 s = —I, 
which allows us to reduce the parameter space for the 
RG flow. We expect that the superconductor-insulator 
quantum phase transition is not sensitive to details of 
interactions and, therefore, concentrate on the simplest 
case of broken spin-rotational symmetry, where we have 
a two-parameter RG flow (for couplings t and 7 c). 

In Fig. [ 19 ^ we plot schematically the expected phase 
diagram for physical electrical resistance p (that can dif¬ 
fer from the NLSM coupling t as discussed in Sec. VI) 


and a parameter 7 c (“generalized superconducting inter¬ 
action” ) characterizing the superconducting correlations 
in the system. On the mean held level, this parameter is 
just equal to 7 c, but beyond the mean-held description it 
also rehects order-parameter huctuations and thus char¬ 
acterizes the overall superconducting coherence (hence - 
tilde), diverging when the true 2D superconductivity is 
established. The superconducting hxed point is then lo¬ 
cated at p = 0 and 7 c = — 00 . As mentioned in Sec. 
jV] above, the SIT is most likely governed by the hxed 
point at p = p* ^ I and 7 c = — 00 , see Fig. [^a). The 
existence of such hxed point is compatible with the RG 
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(a) yJ^l + 7l 



FIG. 19: (Color online) Sketch of possible RG flows for the 
SIT with the unstable fixed point at p = p, and | 7 c| = oo (a) 
and with the metallic phase at 7 c = oo (b, left panel) or the 
fixed point at p ~ 1 and | 7 c| ~ 1 (b, right panel). 


flows shown in Figs. [T]and[^ The SIT occurs through 
the separatrix connecting this fixed point with the trivial 
clean non-interacting fixed point p = 7 c = 0. 

Two more possibilities compatible with Figs. [^and[^ 
are as follows: 


(i) the fixed point at p = oo and 7 c = — oo is unstable 
and the SIT fixed point is located at p ^ I and 
| 7 c| ^ 1, see the right panel in Fig. 19 (b); 


(ii) the SIT fixed point is located at p = oo and 7 c = 


In both these cases the flow towards the superconduct¬ 
ing fixed point would occur in the presence of strong 
superconducting correlations for an arbitrary high resis¬ 
tivity. The existence of the SIT in granular systems (or 
2D Josephson-junction arrays) provides a strong evidence 


against such a scenario. Indeed, in such strongly inho¬ 
mogeneous systems the superconductivity is established 
locally, but sufficiently strong disorder prevents vanishing 
of the total resistance. 


At the same time, our one-loop RG analysis does not 
exclude the possibility of existence of an intermediate 
metallic phase around p ^ 1 at 7 c = —oo, e.g. as shown 
in the left panel of Fig. 19 (b). This metallic phase might 
be governed by either a segment of fixed points or by an 
attractive metallic fixed point. In both cases, this inter¬ 
mediate metallic phase would resemble a so-called “Bose 
metal” mentioned in the Introduction. Experimental ev¬ 
idence for existence of such a phase was reported in liter¬ 
ature. In the case of preserved spin-rotation invariance, 
we do not see any indications for such a scenario. For the 
spin-orbit class, our theory does suggest a critical-metal 
phase with resistivity of the order of resistance quantum, 
somewhat similar to the proposed “Bose metal”. How¬ 
ever, in the phase diagram this phase is separated from 
the superconductor by a narrow insulating region, so that 
the expected sequence of quantum phase transitions is SC 
-I - CM - I. More work is needed to prove or disprove 
the possibility of the intermediate metallic phase both 
within the NLSM formalism and experimentally. 

Returning to the phase diagram in Fig. 19 (a), we em¬ 
phasize that the SIT fixed point at p* ~ 1 and 7 c = —oo 
is reached only at infinite RG scale, corresponding to ex¬ 
actly zero T. This implies that in realistic experiments 
(performed at finite temperature) the flow along the sep¬ 
aratrix might seem as a flow towards an insulator. As a 
result, the critical value of resistance i?c(Tmin) inferred 
from the temperature dependence of the resistivity mea¬ 
sured down to finite T^in might be lower than the true 
critical value i?c = (/i/e^)p*, see Fig. 20 Moreover, 
the ‘critical’ values -Rc(Tmm) extracted in such a manner 
from experimental data may significantly differ for differ¬ 
ent values of Tmin- This example is typical for the two- 
(and more) parameter scaling and demonstrates that the 
“non-universality” of the critical resistance of the SIT 
might be an artifact of the interpretation of the data ob¬ 
tained for finite temperatures. 

It is also worth noting that in the case of preserved 
spin-rotational symmetry (realized in many SIT exper¬ 
iments), the RG flow includes an extra dimension cor¬ 
responding to the triplet coupling 74 . In this situation, 
depending on the initial parameters, the SIT fixed point 
at p* ^ I in the three-dimensional parameter space can 
be reached both from above and from below, thus leading 
to different conclusions on the value of i?c obtained at fi¬ 
nite T„iin in different samples or settings. Moreover, the 
extra dimension in the phase diagram in this case might 
result in a non-monotonic temperature dependence of the 
resistivity at the SIT, as observed in some experiments. 


Let us emphasize that the phase diagram in Fig. [T9|is 
obtained within the framework of the so-called “fermionic 
mechanism” of the SIT (associated to the works by 
Finkel’stein |32I42I47| ). There is a popular misconcep¬ 
tion in literature stating that, contrary to the “bosonic 
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FIG. 20: (Color online) Sketch of a typical dependence R{T) 
near the SIT. The black dashed curve demonstrates how the 
maximal resistance approaches the critical one with decrease 
of T. The solid green curve indicates the separatrix. The 
dashed green line demonstrates approximation to the critical 
resistance from the curves available for T > Tmin (see text). 


mechanism” |33j yielding p* ~ 1, the fermionic mecha¬ 
nism predicts p* <C 1 for a system with Coulomb interac¬ 
tion. Our analysis demonstrates that this is not so: the 
SIT within the “fermionic mechanism” is governed by a 
fixed point with the critical resistivity p* ^ 1. The confu¬ 
sion might arise if one neglects localization effects (due to 
interference and interactions) in the analysis of T,. as in 
Refs. |32l47j . Indeed, for very weak disorder (tg ^ 7co) 
the renormalization of t is negligible, see Figs. [^and|^ 
where the red SC curves appear to be almost horizon¬ 
tal for small to. However, with increasing disorder and 
approaching the separatrix (tg oc y^g <C 1) the red SC 
curves become more and more pushed into the region of 
p ^ 1 at low temperatures, and the critical resistance Rc 
is of the order of Rq. In other words, temperature depen¬ 
dence of resistivity for SC curves close to the transition 
reaches at the maximum a value ~ Rq before dropping 
down with further lowering of the temperature. This is 
very well illustrated by the resistivity plots in the present 
paper. We thus reiterate that both fermionic and bosonic 
mechanisms predict the SIT governed by a fixed point at 
p* ^ 1. 

VIII. SUMMARY AND CONCLUSIONS 

To summarize, in this paper we have explored by 
means of the RG approach the interplay of superconduc¬ 
tivity, interaction, and localization in 2D quantum sys¬ 
tems. The focus has been put on the SIT in thin films. 
Our main results are as follows. 


have derived the full set of one-loop (in disorder 
strength t <C 1) RG equations for a 2D disordered 
interacting system, Eqs. (23) - (26). Formally, these 
RG equations are valid for arbitrary interaction 
couplings, including an arbitrary strong amplitude 
7 c in the Cooper channel. The range of the applica¬ 
bility of one-loop RG equations has been identified 
as the domain t max{l, lyd} 1: beyond this range 
higher loops become important when approaching 
to the insulator t = oo, oi else, to the supercon¬ 
ducting instability, 7 c = —oo. 


2. We have employed the RG framework to explore 
the structure of the phase diagram at zero mag¬ 
netic field. The analysis of RG equations has been 
performed for systems both with and without spin- 
orbit interaction (see Figs. 0 HI n 0 § [m and 


13). Furthermore, the cases of short-ranged and 


long-ranged Coulomb interaction have been inves¬ 
tigated. In general, the phase diagram of a 2D dis¬ 
ordered interacting system is determined by multi¬ 
parameter scaling and depends on the symmetry of 
the problem. 


3. The enhaneement of 2D superconductivity by An¬ 
derson localization m has been confirmed for the 
short-range case. We have identified the parame¬ 
ter regions where the superconductivity is enhanced 
by localization, both for the cases of preserved and 
broken spin-rotational symmetry (see Figs. 
andfl^. 


4. In the case of preserved spin-rotational symme¬ 
try, the RG flow describes the three-parameter (for 
Coulomb repulsion) or four-parameter (for short- 
range repulsion) scaling, thus rendering the phase 
diagram multidimensional, with nontrivial fixed 
points appearing. In particular, a ferromagnetic 
phase develops with the metallic temperature be¬ 
havior of the resistivity in a range of temperatures 
above the ferromagnetic instability (see Fig. |^. 
This behaviour of resistivity in the ferromagnetic 
phase may be confused with a tendency to super¬ 
conductivity in experiments. 


5. The presence of spin-orbit coupling (which removes 
the triplet interaction channel and converts weak 
localization into antilocalization) strongly affects 
the overall phase diagram. Two of this changes 
are fully expected. First, the spin-orbit coupling 
eliminates the ferromagnetic phase. Second, in the 
case of short-range interaction, a super metal phase 
emerges. What is much more intriguing, our results 
indicate that a critical-metal phase with resistivity 
of the order of resistance quantum Rq may arise 
(see Figs. and [^. This phase bears certain 
similarity with a “Bose metal” phase, evidence for 
which has been found in some experiments. 


1. Within the non-linear sigma model formalism, we 


6 . We have evaluated the temperature dependence of 
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the electrical resistivity, p(T), for given bare (high- 
temperature) couplings down to the lowest temper¬ 
atures of the applicability of the one-loop RG ap¬ 
proach. In this temperature range the resistivity 
is dominated by the NLSM coupling t taken at the 
length scale Lt- At lower temperatures, in the close 
vicinity of T^, the electric resistivity is controlled by 
contributions due to inelastic processes (which are 
not automatically included by the RG procedure 
but rather require an additional calculation once 
the RG has been stopped by temperature). These 
corrections lead to a strong suppression of p{T) (see 
Fig. 14). 


7. We have stud ied the magnetoresistance near the 
SIT (Sec. |VI[). Both orbital and Zeeman effects of 
the magnetic field have been included in the uni¬ 
fying RG scheme complemented by the analysis 
of fluctuation corrections near the superconduct¬ 
ing transition. A non-monotonous magnetoresis¬ 
tance, see Fig. [T5|has been predicted, with a max¬ 
imum near the critical field H(,, in agreement with 
experimental observations. The magnetoresistance 
becomes progressively stronger with lowering tem¬ 
perature and becomes giant as T —0, as also seen 
in experiments. 


8 . We have further discussed in Sec. Eni the lim¬ 
itations and generalizations of our approach as 
well as comparison between our theory and exper¬ 
iments. In particular, we have analyzed in detail 
how the screening of Coulomb interaction by sub¬ 
strates with high dielectric constant and by exter¬ 
nal metallic gates can be taken into account within 
our framework (see Fig. 17). 


9. This consideration allowed us to propose specific 
sample designs for experimental observation of the 
superconductivity enhancement. Also, we argued 
that our results for enhancement of the super¬ 
conductivity are in qualitative agreement with 
the experimental observations in layered material 
Li,,ZrNCl. 


10. We have discussed the possible overall structure of 
a generic SIT phase diagram (see Fig. 19) and 
the implications of our findings for the experimen¬ 
tal verification of the universality of this quantum 
phase transition (see Fig. 20). In particular, we 
have shown that, in contrast to a popular belief, 
both “fermionic” and “bosonic” mechanisms of the 
SIT have to do with the fixed point characterized 


by R' 


Rq. 


Our findings are in a qualitative agreement with most 
of the experiments on 2D superconductivity in disordered 
films. In particular, our approach explains a seemingly 
non-universal behavior of the critical resistance Rc found 
experimentally in different systems. Further, the analysis 
of the RG equations in a magnetic field is in accord with 


the experimentally observed non-monotonic behavior of 
the magnetoresistance near the SIT. 

The detailed analysis of the temperature dependence 
of the resistivity in the close vicinity of the “classical” 
(flnite-T) superconducting phase transition governed by 
the Berezinskii-Kosterlitz-Thouless physics is a subject 
of future work [^. The effect of disorder and interac¬ 
tion on the topological sector of the theory (vortices in 
the order-parameter phase) is expected to be increasingly 
more pronounced upon approaching the SIT. A related 
line of future research is devoted to the study of the 
zero-T resistance of a disordered superconducting film as 
a function of the sample size. Further, the RG approach 
developed here will be employed for studying the local 
tunneling density of states (including its mesoscopic fluc¬ 
tuations) in the presence of superconducting correlations 
[45] as well as the fluctuations of the superconducting or¬ 
der parameter near the SIT, as measured in experiments. 
Finally, the role of strong local superconducting fluctu¬ 
ations in physics on the insulating side of the SIT still 
requires further study. 
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Appendix A: Background field renormalization of 
the nonlinear sigma model action 


In this appendix we present details of the derivation 
of one-loop renormalization of the NLSM action 0 with 
the help of the background field renormalization. Let 
us separate the matrix field Q into the “fast” (Q) and 
“slow” (Qo = Tq^KTq) modes as 

Q ^ To-iQTo. (Al) 

We assume that the “fast” matrix field Qnm has non¬ 
trivial structure in the Matsubara space for frequencies 
|n|, |m| < Amax whereas the “slow” matrix field Tq has 
non-trivial structure only for smaller frequencies: 


(^o)nm 


n; 1^1 ^ ^niax and |?7l| ^ Umax? 

\n\ ^ Umax or |m| < Umax- 

(A2) 

In what follows we assume that A„iax ^^rnax ^ 1- 
It is convenient to rewrite the sum of interacting terms 
S[^, S[:} and in NLSM Q as 


"int 




drTrj Tr 


7 “ (1 


tr 


7 “ O 

^n.ri ^ 


an r,j 


(A3) 
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where 


J, 


n,rj 


and 


InUj, r = 0,3, j = 0,1,2,3, 
LnUj, r = 1,2, j = 0,1,2,3. 


(A4) 


For the perturbative (in l/g) treatment of the corre¬ 
lations of the fast fields in Eq. (All) we shall use the 


square-root parametrization of the fast fields 


Q = W + K\/l - 1F2, 


W = 


0 w 
w 0 


(A12) 


r,o = (-l)’T., r,, = -(-l)’Tt, r = 0,3, j = 1,2,3, 
Tro = Tc, = 0, r = 1, 2, j = 1, 2, 3. (A5) 

The effective action for the slow Qq field is given by 


5eff[Qo] = In y V[Q] exp S[T-^QTo] (A6) 

with 

5'[7o ^QTq\ = *S'[Qo] + 'S'[Q] + Ocr + Oint + Or), (A7) 
where 


o, = oi« + o?» + o?>'2 

0... = 0<”■‘+oL■>■* + O^S■'+oS'^ (A8) 

Here we introduce the following terms {5Q = Q — K) 


OW =-I JdrTrA6Q\/6Q 

0^2).1 fdrTiASQAA 


OS")’" = -^ JdrTiASQASQ 
ttT 


0Lt' = -^EE drT^.Tr tr 


ocn r,j 


ttT 


oLE = -AfEEpr.,Tr 


j:,rjSQ 


Ocn r,j 


ttT 


oLr = - EEEKr ., T , 


Ocn r,j 


ttT 




l ( 2).2 __ 

4 

Ocn r,j " 

On = AttTZc^ drti'AnSQ, 


Jn,rjQo 


K,rj6Q 


tr 


tr 


tr 




K.rjSQ 


Al^^SQ 


(A9) 


where 


We adopt the following notations: = Wnin 2 nnd 

bF„ 2 ni = with TT-i ^ 0 and 712 < 0. The blocks w 

and w (in Matsubara space) obey 

vj = -CuFc, w = -Cw*C. (A13) 


The second equality here implies that in the expansion 
K^n 2 = Y.r 3 iKfu 2 )rjtrj souiB of the elements (w“f„Jrj 
are real and some are purely imaginary. 

Expanding the NLSM action ([^ to the second order in 
IF, we find the following propagators for diffusive modes. 
The propagators of diffusons (r = 0, 3 and j = 0,1, 2, 3) 
read 


[Wrj{q)]nlnl[Wrj{-q)]t^Z) = 

327rrr, 


X'Dq[ini2) 


, (A14) 


where Hf 2 = ^ni ~ £n 2 = 27rTrzi2 = 2TrT{ni — 712 ). The 
standard diffuson propagator is given as 

X>“^(7a;„) = + 16Z,^|a;„|/5. (A15) 


The diffusons renormalized by interaction in the singlet 
{T)q°\uj) = 'D^{lu)) and triplet = V^q\uj) = 

'Dq^\uj) = I?*(a;)) particle-hole channels are as follows 


[I?^(fa;„)]-i =q^ + + r«)|a;„|/g, 

[V\{iujn)]-^ = 9 " + 16(^c. + rt)|w„|/5. (A16) 

The propagators of singlet cooperon modes (r = 1,2 and 
j = 0) can be written as 


647rTz 


’^147^32 


xCq(i^\2) 


9 




(A17) 


A = ToVT^\ A, = ro[7;,To-i], Tp-i]. 

(AlO) 

In the one-loop approximation the effective action 
<S'eff[Qo] can be obtained by expansion of ^'[rQ'^QTo] to 
the second order in A^ and Then, we find 

^eff[Qo] = S[Qo] + (O.) + (Oint) + (Or,) 

d--^{{Oa + Oint + On)^), (All) 

where the average (...) is taken with respect to ac¬ 
tion Q. 


where £"12 = Em + £7121 Cqi'icOn) = 'Dqii^n)-, and the fluc¬ 
tuation propagator ( 7 c = Tn/Z^n) 


Cq ^{iujn) = 7c ^ + In ^ - V' 


Dq^ + |tLi„| + A' 1 
FdT 2 


+ V' 



(A18) 


Here D = q/il&Zm) stands for the diffusion coefficient, 
4>{z) denotes the di-gamma function, and A = 

(A' = 47rTr7inax) determines the ultra-violet (infra-red) 
for the fast modes. The propagators of triplet cooperons 
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(r = 1, 2 and j = 1, 2,3) are insensitive to the interaction 
and are as follows 


side of Eq. (All |. We start from corrections to the con¬ 
ductance: 




nina 

^ ^712-714 (A19) 


In general, each term in the right hand side of Eq. 
(All) produce contributions which cannot be expressed 


in terms of Qo only. However, all such contributions can¬ 
cel in the total expression (All). Therefore, we will not 


list them below. Expanding SQ in series of W according 
to Eq. (A12) and performing averaging with the help of 
Eqs. (A14) and (A17), we obtain the action 0 for the 
slow field Qo but with 

5(A) 5(A') = g{A) + Sg, 

rs,t,c(A) —>■ rs_t_c(A^) = rs^t_c(A) -i- 

Z^A) ^ Z^iA') = Z^{A) + (A20) 

The contributions 5g^5Ts^t,c^ and SZ^^ take into account 
the effect of integration over the fast modes from the 
ultra-violet energy scale A down to new ultra-violet scale 
A'. 

Below we list the different non-zero contributions to 
5g,5T and 5Z^jj from each term in the right hand 


J q,(jJn ^ 

3 

+ (A21) 

j=0 


(0(2),2) ^ 5g(2).2 ^ _4 ^c,(0). (A22) 

J a 


Here we use the following notation. 


and 




xQ(^Dg^ + \un \ — A^), 


(A23) 


and 


(fq 


0 (A - Og2) 0 I^Dq^ - A'), (A24) 


J,~ J ( 27 r )2 

where 0(a:) stands for the Heaviside step function. Next, 


O 


(2).2 


2 L 


O. 


(1),2 


^ = -2 


^-/j'Dq{iuJn)'D[^'^{iU}n) [l “ 2q^T>q{iuJn)] 




-4 f Cq{iUJn)Cl{iUJn) [l “ 2q^Cq{iWn)\ ■ 

J q,UJn 


In total, we find 


5(A')=5(A)-4 /c,(0)+4 / g2 


2Cq{iUJn)Cq{iuJn) + ^ ^j'Dq{iuJn)'D^q^\iUJ„) 
j=0 


The renormalization of Z^^ is described by the following terms: 


(A25) 


(A26) 


(O^) ^ / \2Z^C^^{iuJn)Cq{iuJn) + ^ ^qVq{^COr^)Vlf'>{i0Jn) 


1=0 


(A27) 


O 


(2).2 


2 L 


O 


( 1 ). 2 ' 


= _ / J ^ r, - l]dnDqiiuJn) + 2Z^Cq{iuJn)dnCq{lUJn) L 

(A28) 


where i9„/(za;„) = ((//16)9//9a;„. Combining Eqs. (A27) and (A28) together, we obtain 


Zu]{A') — Ztjj{A) -\ -(rs-|-3r4) f 'Dq{iuJn) - Zi^ f Cq{iUJn)dnCq{iuJn) — Cq{iUJn)dnCq{iuJn) ■ (A29) 


-qyvs^fij I I '^qy 

9 J q,UJri 9 Jq,U)ri 


The corrections to the interaction amplitude Tg are as follows: 


(oLt') ^ srfj’l = T, / 2Z^C2(za;„)/:,(za;„) + J] (za;„) 

9 J q,UJri ^ 


(A30) 


1=0 






























28 


+ 2 


O, 


(1).2 


n2 


s,int 


:- -(r. + sro [v,{o) 

9 J q 

- [ Cg{0)Cg{0) - j Cl{i0Jn)Cl{lUJn), 

9 Jq 9 J q^ujn 


^(i).1q(i). 2\ ^p(i),i;(i).2 _ _4 r 

'^int ^int / ^ s,int ~ 

' 9 J (},^n 


2Zu:Cq{iuJn)Cq{iuJn) + '^T jVg{iu;n)'D^q^\iuJn) 
j=0 


In total, we find 


r,(A') = r,(A) - -(r, + 3rt) [ Vg{Q) - -Z^ /'c,(0)/:,(0) - -Z^ f C^(iu;„)£^^(iu;„). 
9 Jq 9 Jq 9 J q.uJn 


f q J q 

The corrections to the interaction amplitude Ft can be listed as follows: 


(oLt') ^ = -tJ 2Z^C^(tu;„)£g(tu;„) + ^ 

9 


o 


(2).2 


2 L 


o 


(1).2 


n2 


^rlfnf = -"(r* -o) /12,(0) + / c,(0)£,(0) - 


I ^qy^JJ^qy'JJ — -1^ t I ‘^^q‘ 

9 Jq 9 Jq 9 Jq,(jJn 


o 


(i).i' 


(l),l;(i).i _ °p2 / 

9 J q^UJn 


^ t/mt ■*- 


V^JiuJn) - VfiiuJn) 


K2f D 


- -r 


9 ‘ V47rT 


£q{i'^n)'4^ {^9^q,i\ujj^\) 


oToT) ^ = )r, f 

' 9 J q.UJn 


2Z^Cq{iuJri)£q{i(^n) - j'£>qiil^n)'£>q\ii^n) + 4^ (iuj^) 

i=0 


Here 9(^q^i\uj^\ = {Dq^ + |w„| + A')/(47rr) + 1/2.In total, we obtain 

rt(A') = r,(A) - ^(r, -r*) [Vqio) + T/ [ vl{iuj^) + /'c,(o)£,(o) 

9 J q 9 J q.oJn 9 J q 


H- £tZijj f £q{iuJn) 

9 J q,0Jri 


2z\A^T 


i’" {^q,Au:„\) 


Finally, the corrections to the Cooper channel interaction amplitude Fc are given as 

3 


')( 2),2 

^int 




1 

+ 2 

'^int 


A2).i ^ 4 

c,int g 


i=o 


n 2 




: - -(r.-aro f Vq{o) 

9 Jq 

X>“^ (iWn {iuJn)£q{iuJn) 


- -TsZ^ 


’ Qi^n 


X>g(za;„) +Cg{iuJn) 



ST 


(i),i;(i),i 

c,int 


T' c^q—oi'^^n 



Cq—Q (iu^TTT,), 


(A31) 

(A32) 

(A33) 

(A34) 

(A35) 

, (A36) 

(A37) 

(A38) 

(A39) 

(A40) 


1 

2 


(A41) 
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^(l).l^(l).2\ .p(l),l;(l),2 _ 4 

"^int ! ^ c,int ~ 


/ ^q—o(^^m) 

1 [ 


J q,(jJn 


2Z^Cg{iU}n)j^q{iU}n) 


j=0 


Here we introduce 

In total, we obtain 

re(A') = r,(A) + r, 


2nT 


0(A - Dq^ - \oj^\)e{Dq^ + |a;„| - A'). 

n>0 


(A42) 


(A43) 


/Lq—oi^iUJri — ^) / o(^^m) 


2Za;Cq{iUJn)Cq{it^n) + 


j=0 


--(r, - ar*) [ i?,(o) - h,z^ j 

9 Jq 9 Jq.ujr, 


2^0 (*w„) + Cg{iUJn) 


(A44) 


We emphasize that Eqs. (A29) and (A33l implies that 


4 f r 47rT i 

Z^(A') + r,(A') = Z^A) + r,(A) + -Z,, / Cqiiojr^Yqiiojn) - —C,(zw„) = ZU^) + r«(A). 

9 Jq,UJ-n ^ U J 

(A45) 

Here we employ the following relation: 9„£q(zw„) = [-D/(47rT)]£^(za;„)V''(A’^,i|i^„|)- Also since A' 47rT we use that 
V'' Yq,i\uj^\) ~ (47rT/D)C,(za;n). 

The one-loop renormalization of the parameters of the NLSM action obtained from the background field procedure 
can be summarized as follows: 


5(A')=5(A)-4 /C,(0)+4 / q^ 

9 q 9 q,^n 


2Cg{iuJn)Cq{iuJn) + Y. 

3=0 


(A46) 


2 f f 

^cj(A^) = Zi^(A) H-(Tg -f aTi) / 'Dq{iuj) H-/ Cq{iuJn) Ci{iU3n) +■ Cq{iUJn) 


(A47) 


rg(A') = rg(A) - ^(Tg + 3r0 [ V^w) -7 a f C'iituJr.) + Cq{^u;n) 

9 J q,UJri 9 J q,UJri 


(A48) 


and 


r*(A') = rt(A)-^ frg-Ti-d^' 

9 \ 


4 / 2T^ 

'D‘^{iuJn) Y - i Y 4:Tt Y I . 

/ 9q^uJn 9 \ J 9q^uj 


C {iuj Y 7z)/Iq(zu^^) 




(A49) 


r,(A') = r,(A) - ^ / C,=o(*Wn) - -(r. - ar*) [ V^Wn) + — [ CYoJn)Cq{iuJn) 

^LJ 9 oJn 9 9 q^ijJn 9 9 q^oj^ 


i2rtr, 

9Zu, 


Vq{iuJnYq{iuJn). 


(A50) 


q^tjJn 


To derive RG equations from Eqs. (A46) - (A50) we choose A' = A-|-dA. Then, provided dA/A corresponds 
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to —2dL/L we obtain RG equations (18) - (22). 


Appendix B: Enhancement of in the case of weak 
short ranged interaction 


1. Orthogonal symmetry class 


Expanding RG Eqs. (23) - (26) with n = 3 to the 
lowest order in interactions 7s,t,c we find 


dt n 

_ =f^ 

dy 



Here within our accuracy we neglect the interaction cor¬ 
rections to t in comparison with weak-localization cor¬ 
rection. The matrix 


/I 3 2 \ 

= 1 -1 -2 (B2) 

V1 -3 0 y 


has the following eigenvalues and eigenvectors: 


A = -4 



A' = 2 



If the 7 ^ te rm is neglected, the solution of the linear 
system (Bl) approaches the eigenvector with A = —4, 
i.e., interaction parameters tends to the BCS line — 7 s = 
7 t = 7 c. Le t us expand the vector formed by 7s,t,c in 
eigenvectors (B3): 



Transforming the set of equations (Bl) to the new vari¬ 
ables a, 6, and c, we get 


dt 

dy 

da 

dy 

db 

dy 

dc 

dy 


2 o 

2ta — -{a — b + 2cY, 

O 

—tb, 

2 

—tc — -{a — b + 2c)^. 
o 


(B5) 


Erom Eqs. (B5) it is easy to find 

tiy) = 


to 

1 - ^oy’ 


to 


b{y) = boY-toy) = boj- (B7) 


Since b decreases upon RG flow, it is not important and 
we neglect it in the future analysis. 

Equations for the remaining two variables, a and c are 
coupled. If the quadratic term is neglected, then a in¬ 
creases and c decreases. This suggests that c can be 
neglected. This is confirmed by a more careful analysis 
which shows that, although on the very last interval of 
RG “time” y the variable c starts to increase and becomes 
of the same order as a (i.e. of order unity), this weakly 
affect the RG scale at which this happens (i.e. the tem¬ 
perature of the superconducting transition). Thus, we 
neglect c in what follows. 

We can now easily solve the remaining equation for a. 
We assume the starting value oq to be negative (which 
means that there is attraction in the Gooper channel that 
is supposed to lead to the superconductivity), a = — |a|. 
This is in particular the case when 7c,o is the dominant 
coupling and 7c,o < 0. Then the equation reads 


d|a| 

dy 


2t|Q:| 4“ . 

O 


Solving this equation, we obtain 


a{y) 


2to_2^\ ^ 

|ao|t^ 3t^ 3t y 


(B8) 


(B9) 


Let us analyze the result obtained. Let us first assume 
that |ao| ^ to- Then the second term in brackets in the 
right hand side of (B9) is small compared to the first one 
and can be neglected. 


a-Yy) 


1 

t 




(BIO) 


With increasing RG scale y the resistance t increases to¬ 
gether with the interaction a. If t reaches first unity, 
we get an insulator; if a ^ 1 happens first, we get a 
superconductor. It is easy to see that the second possi¬ 
bility (superconductivity) is realized if |ao| ^ tg- Then 
at the point of divergence of a we have a resistance 
G ~ 3tg/(2|ao|) <C 1. This occurs at y* ~ T _ i, 
i.e. we can estimate the temperature of superconducting 
transition as Tc ~ exp(—2y*) yielding 


Tc 




(Bll) 


Equations (B5) are supplemented by the following initial 
conditions: t(0) = to, a(0) = ag, b{0) = bo, and c(0) = cg 
where 


/ao\ /-1/6 1/2 l/3\ /js,o\ 
Ug = 1/2 1/2 0 7t.o . 

\coJ V V 3 0 I/Sj \jc,oJ 


(B6) 


Here the factor 2 in the exponent originates from a trans¬ 
lation of the length scale into energy (temperature). Un¬ 
der the above assumption |ag| ^ to the s econd term in 
square brackets in the exponential of (Bll) is just a small 

correction to the first one. _ 

The transition temperature (Bll) is much higher than 
the BGS temperature T®®® = (l/r)e“^'^l^‘=’“l, so that 
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the superconductivity is strongly enhanced by disorder. 
The origin of the enhanced superconductivity is in the 
increase of |a| governed by the eigenvalue A = —4 of 
matrix TZo which yields the eigenvalue — (t/2) x (—4) = 2t 
of the linear part of the system in Eqs. (Bl). This is 


nothing but the anomalous multifractal exponent — A 2 
for this symmetry class [52] • (We have in mind the “weak 
multifractality” in 2D.) Therefore, the (multi)fractality is 
the source of the enhancement of the superconductivity. 


By solving (B5) with 6 = 0 and a given by Eq. (BIOI, 


one finds that although |c| decreases initially, eventually 
with increasing RG scale towards y* it becomes of the or¬ 
der of unity: |c( 2 /*)| ^ 1. Therefore, to determine precise 
value of G one has to solve coupled equations for a and 
c [with (6 = 0 )]. 

If |ao| tgi the resistance reaches unity before the in¬ 
teraction becomes strong, and the system is an insulator. 
Finally, if |ao| ^ toi the disorder is not particularly im¬ 
portant, and the transition temperature is given by usual 
clean BCS In the latter case neglecting 6 and c 

is not parametrically justified and leads to an incorrect 
numerical factor in the exponent. 

For the initial values of interaction parameters used in 
the inset of Fig. j^we find that oq = —0.01. Thus, we ex¬ 
pect the regime with Tc ~ for to |ao| = 0.01. In 

the range of to between |ao| = 0.01 and ^|ao| = 0.1 the 
transition occurs at Tc S> At to ~ a/IooI =0.1 

the superconductor-insulator transition is expected. The 
above crude analysis of RG equations linearized in inter¬ 
actions is in good agreement with numerical solutions of 
full RG equations (23) - (26 1 . 


2. Symplectic symmetry class 

Now we consider the symplectic symmetry class, i.e. 
assume that the spin symmetry is completely broken. 
Expanding Eqs. (23) - (26) with n = 0 to the lowest 


order in interactions 7s,t,c we obtain 


— -iv 

dy 2' dy \lcj 2 " V7. 


1, 


0 

H 


2 • (B 12 ) 


Here we neglect interaction corrections to t in comparison 
with weak anti localization. The matrix 


TZs = 


1 2 
1 0 


has the following eigenvalues and eigenvectors: 


A = -l : 


-1 

1 


A' = 2 : 


1/ ■ 


(B13) 


(B14) 


If the 'Yc term is neglected, the solution of the linear 
system in Eqs. (B12) approaches the eigenvector with 


A = — 1, i.e., interaction amplitudes approach the BCS 
(in the absence of 7 *) line 7 s = — 7 c- As in the orthogonal 


case, we can expand the vector formed by 7 s,c in the 
eigenvectors 


-1 2 

1 1 


(B15) 


Transforming the set of equations (B12) to the new vari¬ 
ables a and c, we find 

dt t^ 

dy 2 ’ 

da t 4, ,0 

=-tc-ha + cf. (B16) 

dy 3 


Equations (B16) are supplemented by the initial condi¬ 


tions: t( 0 ) = to I q( 0 ) = oo and c( 0 ) = cq, where 


ao\ _ 6-1/3 2/3 
Co \ 1/3 1/3 


7s ,0 
7c,o 


From the first of Eqs. ( |B16 1 we find 

to 


t{y) = 


1 + 2/io/2’ 


(B17) 


(B18) 


Equations for two variables, a and c are coupled. If the 
quadratic term is neglected, then a increases and c de¬ 
creases. At the later stage of RG the quadratic terms 
leads to enhancement of c. This suggests that c can b e 
neglected for qualitative analysis of RG equations (BI 61 . 


We thus neglect c and keep only a (fully analogously to 
what we have done in the orthogonal case). The resulting 
equation for a reads 


da 

dy 


t 4 o 

= 7 : 0 -— 7:0 

2 3 


We solve this equation with the result 
1/14 4 


t \aoto 36^ 


3^0 


-1 


(B19) 


(B20) 


Provided Oq < 0, the superconducting instability is pos¬ 
sible. The new, different from standard clean BCS be¬ 
havior emerges under the condition |ao| ^ to- Then the 
condition a ~ 1 yields 


G ~ 2 (|ao|to/3)^^^ < to. 


(B21) 


By solving Eq. ( B16| ) with a given by Eq. (B20l, we find 
that although |c| decreases initially, eventually with in¬ 
creasing RG scale it reaches a: c ~ a at t = t*. Therefore, 
to determine precise value of t* one has to solve coupled 
equations for a an d c. 

Equation (B21| yields the following estimate for the 


transition temperature 


Tc ~ 

T 




exp 


\/\oo\to 


> T 


■BCS 


(B22) 
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The constant C is of the order of unity and depends on 
ratio co/ao- The clean BCS result is restored (up to small 
corrections) at |ao| ^ to- For oq < 0 (or, equivalently, 
7c,0 < 7 s,o/2) there is no transition to the supermetallic 
phase with increase of to- 

As in the orthogonal case, the source of the enhance¬ 
ment of the superconducting temperature is in the first 
term on the right hand side of Eq. (B19). The eigenvalue 


t/2 is the anomalous multifractal exponent —A 2 for the 
symplectic symmetry class. Therefore, also in this case 
the (multi) fractality is the source of the enhancement of 
the superconductivity. This enhancement is less efficient 
than in the orthogonal case for two reasons, because of 
antilocalizing behavior that leads to decrease of t and 
therefore weakening of multifractality. 
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